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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  1976 

SUMMARY  OF  INFORMATION  RELATING  TO  GUST 
LOADS  ON  AIRPLANES 
By  Philip  Donely 

SUMMARY 


Available  information  on  gust  structure,  airplane  reactions,  and 
pertinent  operating  statistics  has  been  examined.  This  paper  attempts 
to  coordinate  this  information  with  reference  to  the  prediction  of  gust 
leads  on  airplanes.  The  material  covered  represents  research  up  to 
October  19^7* 


INTRODUCTION 


The  fact  that  all  airplanes  fly  in  rough  air  at  some  time  poses  a 
number  of  problems  relative  to  safe  flight.  One  of  the  most  important 
of  these  problems  is  that  of  designing  the  airplane  structure  to  with¬ 
stand  the  loads  imposed  by  gusts.  The,  three  principal  phases  of.  the 
gust -load  problem  are:  (l)  the  determination  of  the. gust  structure 
(that  is,  the  size,  shape,  intensity,  and  frequency  of  occurrence) , 

(2)  the  reaction  of  any  airplane  to  gusts  of  known  structure,  and 

(3)  the  determination  of  the  operating  statistics. 

No  order  of  importance  can  be  given  to  the  three  phases  of  the 
problem  since  the  final  loads  are  a  function  of  the  gust,  the  airplane, 
and  hew  the  airplane  is  flown.  The  characteristics  of  gusts,  where  they 
occur,  and  the  variation  in  their  characteristics  are  of  fundamental 
importance  because  the  gust  is  the  source  of  the  problem.  Analytical 
end  experimental  work  on  what  happens  to  an  airplane  when  it  strikes  a 
known  gust  is  of  importance  ..since  a  knowledge  of  airplane  reactions 
permits  the  load  calculation  for  any  airplane.  Finally,  operating 
statistics  are  of  importance  in  setting  the  level  of  loading  for  operating 
airplanes . since  they  define  the  gusts  encountered  under  actual  operating 
conditions  and  the  speeds  at  which  the  gusts  are  encountered. 

Although  research  on  gust  leads  has  been  carried  on  for  many  years 
and  many  of  the  results  have  been  incorporated  in  design  rules  (refer¬ 
ence  l),  these  results  have  either  been  issued  piecemeal  or  not  published 
at  all.  Under  these  conditions,  a  compilation  of  available  information 
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within  a  logical  framework  is  needed  and,  for  this  purpose,  the  present 
paper  has  teen  prepared.  The  scope  includes  all  NACA  material  available 
up  to  about  October  1947  arranged  according  to  the  three  principal  phases 
previously  mentioned. 

Because  airplanes  are  used  to  obtain  gust  data  and  the  various 
phases  of  research  are  interdependent,  the  first  section  of  the  paper  is 
a  presentation  and  discussion  of  the  basic  equations  of  the  reactions  of 
an  airplane  as  used  in  gust-load  calculations.  This  first  section  is 
followed  by  sections  covering  gust  characteristics,  airplane  reactions, 
and  operating  statistics  and  by  a  short  concluding  section  on  the  coordi¬ 
nation  of  information  and  the  calculation  of  applied  loads.  The. material 
contained  in  this  paper  is  based  mainly  on  past  work  of  the  NACA,  much 
of  which  has  been  accomplished  through  the  cooperation  of  the  agencies 
and  airlines  listed  in  appendix  A. 


BACKGROUND  AND  BASIC  EQUATIONS  OF 
GUST  -  LOADS  RESEARCH 


In  the  earliest  stages  of  research  on  gust  loads,  an  immediate 
practical  answer  to  the  question  of  appropriate  design  loads  was  needed. 

The  problem  was  approached  by  recourse  to  the  obvious  and  necessary 
expedient  of  measuring  the  accelerations  under  actual  operating  condi¬ 
tions  so  that  some  data  could  be  made  available  to  designers  pending  a 
more  complete  understanding  of  gust-load  phenomena. 

Accelerometers  were  used  in  this  work  although  it  was  evident  that 
accelerations  measured  on  one  type  of  airplane  would  not  be  applicable 
to  other  types,  or  even  to  the  same  type  operated  at  a  different  speed. 

The  main  effects  of  airplane  characteristics  and  airspeed  were  taken 
into  consideration  in  overcoming  this  difficulty  by  utilizing  a  formula 
based  on  the  most  elementary  concepts  of  the  nature  and  action  of  gusts. 
This  formula,  called  the  "sharp-edge-gust  formula,"  was  applied  In  such 
a  way  that  the  accelerometer  data  were  reduced  to  effective  gust 
velocities,  the  term  effective  being  employed  to  indicate  the  fictitious 
nature  of  these  gust  velocities. 

Since  that  time,  several  advances  have  been  made  so  that  at  present 
three  concepts  of  evaluating  airplane  reactions  are  being  utilized: 

(l)  the  very  simple  sharp-edge -gust  formula,  (2)  extended  equations 
taking  into  account  the  gust  shape  and  additional  airplane  characteristics, 
and,  finally,  (3)  the  sharp -edge -gust  formula  modified  by  an  alleviation 
factor.  These  relations  and  concepts  are  presented  subsequently  and  form 
the  basis  of  NACA  research  on  gust  loads. 
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SHARP  -EDGE  -GUST  FORMULA 


The  simple  sharp -edge -gust  formula  was  derived  on  the  basis  of 
numerous  simplifying  assumptions  which  include  the  following: 

1.  The  gust  is  sharp-edged  in  the  direction  of  flight  and  represents 
an  instantaneous  change  in  wind  direction  or  speed. 

2.  The  gust  velocity  is  uniform  across  the  span  of  the  airplane  at 
any  instant  of  time  or  position  of  the  airplane  in  space. 

3-  The  gust  direction  is  normal  to  the  lateral  axis  of  the  airplane. 

if.  The  airplane  is  in  steady  level  flight  prior  to  entry  into  the 
gust,  and  the  airplane  flight  path,  attitude,  and  ground  speed  are  not 
affected  by  the  action  of  the  gust  on  the  airplane.  (This  assumption 
might  be  visualized  by  imagining  the  airplane  to  be  driven  along  a  fixed 
track  through  the  gust  so  that  the  airplane  has  no  resulting  motions  due 
to  the  action* of  the  gust.) 

5.  The  primary  effect  of  encountering  a  gust  is  to  change  the  lift 
on  the  airplane. 

6.  The  lift  increment  of  the  horizontal  tail  due  to  the  action  of 
the  gust  is  negligible  as  compared  to  the  wing  lift  increment. 

7-  The  lift  coefficient  of  a  wing  is  a  unique  function  of  angle  of 
attack  and  independent  of  time. 


For  the  assumptions  noted,  the  expression  for  lift  qr  load  factor 
for' a  gust  of  any  angle  is:  _  vj  j- 
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(See  appendix  B  for  definitions  of  symbols.)  The  term  in  brackets 
represents  the  effect  of  the  increased  resultant  speed  and  the  term  in 
parentheses, .  the  lift  in  steady  level  flight  plus  the  lift  increment 
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due  to  the  change  in  angle  of  attack.  If  the  particular  gust  angles 
of  90°,  0°,  and  l80°  are  considered,  then  the  expression  reduces  to 


for  P  =  90°  and 


|  =  l+  An  =  l  +  ~!  ^IJV 
W  da  W  2 
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for  P  =  0°  and  l8o°. 

The  so-called  sharp-edge-gust  formula  is  the  equation  for  P  =  90°  for 
sea-level  density  and  is  generally  written 


/'  dCL 
P°Ue7e  dct 


An  = 


W 

2s 


The  effect  of  gust  direction  on  the  acceleration  increment  is 
shown  in  figure  1  for  airplanes  with  wing  loadings  from  15  to  60  pounds 
per  square  foot  flying  at  200  miles  per  hour  through  a  gust  velocity 
of  15  feet  per  second.  The  ratio  of  the  acceleration  increment  for 
given  values  of  wing  loading  to  the  maximum  value  for  each  wing  loading 
is  plotted  against  the  gust  angle  P-  Figure  1  indicates,  as  do  refer¬ 
ences  2  and  3>  that  the  acceleration  is  a  maximum  for  angles  close  to  90u 
figure  1  also  indicates  that  the  vertical  or  near-vertical  gust  is  k 
to  15  times  as  effective  in  producing  acceleration  as  the  horizontal 
gust.  Therefore,  in  reducing  acceleration  data,  the  assumption  is 
usually  made  that  the  significant  accelerations  caused  "by  gusts  result 
from  vertical  gusts,  that  is,  P  =  90°. 

The  sharp-edge-gust  formula,  sometimes  used  with  a  correction 
factor,  requires  fairly  rigid  definition  of  the  quantities  to  he  substi¬ 
tuted.  The  equation  is  used  for  general  research  studies  where  masses 
of  acceleration  data  are  to  he  reduced  and  compared  for  evaluation  of 
past  airplane  gust-load  experience  and  for  the  prediction  of  load  exper¬ 
ience.  For  research  studies,  the  actual  weight  of  the  airplane  is  used 
to  determine  the  wing  loading  since  the  data  from  different  sources 
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are  to  be  compared.  For  load-evaluation  studies  of  operating  conditions, 
the  gross  weight  is  used  because  it  yields  conservative  values  of  the 
effective  gust  velocity.  In  all  cases,  the  acceleration  increment  is 
defined  as  the  acceleration  minus  one.  The  other  quantities  have  been 
standardized  as  the  result  of  instrumental  limitations  so  that  sea-level 
density  is  used  with  indicated  airspeed  (assumed  to  be  equal  to  V  )  and 
gross  wing  area.  The  slope  of  the  lift  curve  is  generally  obtained  from 

the  formula  ,'k  kV  t 

1  ^  .• 

jA  /  ^  .v-’ 


airplane 

For  preliminary  evaluation  of  acceleration  data,  sometimes  a  lift-curve 
slope  must  be  assumed  since  the  wing  configuration  is  not  known.  In 
these  cases,  a  value  of  4-5  per  radian  has  usually  been  selected. 

In  the  case  where  loads  are  to  be  calculated  by  utilizing  effective 
gust  velocities  in  the  sharp-edge -gust  formula,  the  quantities  used  must 
be  consistent  with  those  described.  If  different  procedures  are  followed 
the  answers  obtained  obviously  are  incorrect  in  proportion  to  the  devia¬ 
tion.  The  simple  formula  with  its  simplifying  assumptions  limits  the 
operations  that  can  be  performed  on  derived  data. 

EXTENDED  GUST  EQUATIONS 


(?) 


yCX 

6k'  '  5 


wing 


A  +  2 

\ 

\ 


The  restrictive  assumptions  of  the  sharp-edge-gust  formula  caused 
little  concern  until  it  was  to  be  used  for  design  calculations  of  gliders 
and  airplanes  whose  wing  loading  and  other  characteristics  differed 
widely  from  those  of  the  airplanes  that  were  initially  used  to  establish 
the  effective  design  gust  velocities.  The  assumptions  of  infinite  gust 
gradient,  steady  lift,  and  no  vertical  motion  had  to  be  eliminated  in 
the  derivation  of  the  equations  for  airplane  response  to  a  gust.  The 
basis  for  the  final  solutions  was  the  work  of  Kiissner  reported  in  refer¬ 
ence  2,  which  resulted  in  the  analytical  study  presented  in  reference  4. 

For  the  application  of  Kussner's  work  to  the  problem  of  gust  loads, 
the  assumptions  already  cited  in  the  derivation  of  the  sharp-edge-gust 
formula  have  been  continued  with  the  following  exceptions:  (l)  The  gust 
velocity  has  been  assumed  to  be  uniform  across  the  span  of  the  airplane 
at  any  instant  and  to  increase  linearly  with  distance  in  the  direction 
of  flight  until  the  maximum  gust  velocity  U  is  attained^  (2)  the  air¬ 
plane  can  rise  but  does  not  pitch  under  the  action  of  the  gust  and  is  in 
steady  level  flight j  (3)  the  unsteady-lift  functions  of  Kussner  and 
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Vagner  (reference  2) ,  which  were  derived  for  the  two -dimensional  wing, 
are  assumed  to  he  made  applicable  to  the  finite  wing  by  substituting  the 
slope  of  the  lift  curve  for  the  finite  wing  in  place  of  the  slope  of  the 
lift  curve  for  the  infinite-aspect -ratio  or  two-dimensional  wing;  and 
(4)  except  for  the  shorter  gust-gradient  distances,  the  acceleration 
peat  was  assumed  to  coincide  in  space  with  the  position  of  the  point  of 
first  attainment  of  maximum  gust  velocity. 

The  symbols  used  in  the  equation  for  the  gradient  gust  are  shown 
in  figure  2,  and  the  unsteady-lift  functions  used  are  shown  in  figure  3, 
taken  from  reference  5*  It  is  indicated  that,  while  a  sudden  angle-of- 
attack  change,  such  as  is  obtained  by  rotating  or  pitching  a  wing,  yields 
one  variation  of  lift  ,  the  gradual  immersion  of  the  wing  into  a  gust 

yields  a  different  curve  CT  • 

Lg 

The  equation  of  vertical  motion  can  be  writtep  as 


t 


<4  ^  q  (tl  - 1) 

2  da  Lg  x 


u 

V 


dt 


where  the  first  integral  is  the  force  due  to  the  gust  and  the  second 

Integral  results  from  the  vertical  motion  of  the  airplane.  Since  the 

unsteady-lift  functions  Ct  and  CA  depend  on  chord  lengths  from  the 

a  ^g 

start  of  a  disturbance,  the  equation  is  generally  transformed  so  that  the- 

Y 

distance  traveled  s  =  —  t  rather  than  t  is  the  independent  variable. 
The  equation  becomes 


CT  (Sl  - 


fSl 

s) ds  -  —  I 
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from  which  the  acceleration  ratio 


is  seen  to  he  a  function  of 


the  rate  of  development  of  transient  lift,  the  mass  parameter,  and  the 
shape  of  the  acceleration  curve. 


The  preceding  equation  has  been  solved  by  Fredholm's  method  (refer¬ 
ence  6)  and  the  solutions  obtained  are  of  the  form 


where  C  and  B  are  constants  representing  solution  of  the  integrals 
in  the  equation  and  whose  values  depend  only  on  the  distance  penetrated 
into  the  gust.  Figure  4  gives  the  solution  for  a  range  of  gradient 
distance  and  mass  parameter.  Since  the  unsteady-lift  functions  Ct 

g 

and  CL  given  in  figure  3  were  approximated  for  distances  above 

about  8  chords,  the  curves  for  the  larger  gradient  distances  in  figure  4 
are  also  approximate.  The  dash  line  in  the  figure  represents  a  rough 
limit  of  applicability  of  the  gradient -gust  solution  obtained  by  this 
method.  Above  the  dash  line,  the  solutions  for  finite  gradients,  are 
not  valid,  and  the  curve  for  H  =  0  is  to  be  used.  The  dash  line  was 
determined  on  the  assumption  that  no  gradient  gust  has  a  peak  accelera¬ 
tion  occurring  earlier  than  a  gust  with  zero  gradient  distance- 

Figure  4,  together  with  the  sharp-ed^e-gust  formula,  can  be  used 
to  evaluate  flight  records  to  obtain  the  true  gust  intensity.  The 
gradient  distance  H  =  Ytj_  and  airplane  mass  parameter  pig  permit  the 

determination  of  the  acceleration  ratio  from  figure  4.  The  substitution 
of  pertinent  values  in  the  sharp-edge -gust  formula  yields  an  effective 
gust  velocity  which,  when  divided  by  the  acceleration  ratio,  gives  the 
true  gust  velocity.  In  this  case,  true  airspeed  and  the  actual  air 
density  are  to  be  used  in  the  sharp-edge -gust  equation. 

In  the  evaluation  of  time-history  data  of  airspeed,  altitude,  and 
acceleration  to  obtain  gust-gradient  distance  and  true  gust  intensity, 
the  following  methods  have  been  used:  The  acceleration  increment  is. 
taken  as  the  acceleration  increment  at  the  peak  minus  the  steady  value 
prior  to  the  start  of  the  gust.  In  the  determination  of  the  gradient 
distance,  the  time  from  zero  to  peak  acceleration  and  the  true  airspeed 
is  used.  In  expressing  the  gradient  distance  in  chords,  the  mean 
geometric  chord  is  generally  used. 
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GOST  ALLEVIATION  FACTOR 


Much  of  the  data  obtained  from  commercial  operations  with  the  NACA 
V-G  recorder  and  some  gust  data  from  special  investigations  are  not 
amenable  to  analysis  by  the  extended  equations  of  the  previous  section 
because  of  limited  information  or  insufficiently  defined  acceleration 
peaks,  and  so  means  to  account  for  the  main  effects  of  alleviation  had 
to  be  established.  Since,  for  much  of  the  work,  the  level  of  gust 
intensity  has  been  set  by  the  airline  data  and  special  investigations, 
the  results  of  the  extended  equations  have  been  used  mainly  to  adjust 
measurements  of  effective  gust  velocity  for  the  effects  of  unsteady  lift 
and  vertical  motion.  The  gust  alleviation  or  correction  factors  so 
obtained  were  consequently  based  on  a  knowledge  of  the  average  gust 
characteristics  and  the  response  of  the  conventional  airplanes  on  which 
considerable  data  had  been  accumulated.  • 

The  gust  alleviation  factor  is  defined  as  the  relative  response  of 
two  airplanes  encountering  the  same  gust  with  the  gradient  distance 
defined  in  chords.  If  a  representative  gust  size  is  determined  and  a 
particular  airplane  selected  as  being  representative  of  the  airplanes 
used  for  gust  measurements,  then  the  effective  gust  velocities  are 
related  by  the  acceleration  ratios.  Thus,  if  all  gust  measurements  were 
made  on  airplane  1  and  the  size  of  the  average  gust  were  known,  then 
effective  gust  velocities  from  the  accelerations  experienced  by 
airplane  2  are  related  to  those  for  airplane  1  by  the  expression 


The  effective  gust  velocity  representing  the  conditions  is  Ue  or,  if 
measured  with  airplane  2,  is  ^ 
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The  ratio  of  acceleration  ratios  is  the  gust  alleviation  factor  and  may 
he  determined  hy  calculation  or  ty  experiment. 

For  load  calculations,  a  specific  alleviation  factor  K  was  derived 
and  is  shown  in  figure  as  reproduced  from,  reference  1.  The  factor  K 
is  calculated  with  the  Boeing  B-247  airplane  as  a  reference  and  the 
assumptions  were  made  that  wing  loading  is  proportional  to  mass  parameter 
and  that  the  effect  of  pitch  on  the  gust  load  increment  is  the  same  for 
all  airplanes.  The  curve  is  calculated  as  the  acceleration  ratio  for  any 

/w 

airplane  to  the  acceleration  ratio  for  the  Boeing  B-2V7  airplane  (-  =  1 6, 
\  ■  \s 
C  =11/  when  both  airplanes  traverse  a  gust  with  a  gradient  distance 
of  10  chords.  The  resulting  curve  is  based' on  the  assumption  that  all 
airplanes  are  similar  to  the  Boeing  B-2^7  airplane  but  have  different 
wing  loadings. 


THE  STRUCTURE  OF  ATMOSPHERIC  GUSTS 


The  determination  of  the  characteristics  of  atmospheric  gusts 
pertinent  to  the  airplane,  their  variations,  and  their  frequency  of 
occurrence  in  the  atmosphere  is  of  fundamental  importance.  Considerable 
research  has  been  undertaken,  therefore,  to  cover  various  aspects  of  the 
gust -structure  problem.  The  research  performed  includes  flight  tests  in 
the  neighborhood  of  Langley  Air  Force  Base  over  a  period  of  about  10  years 
and  special  investigations  made  in  conjunction  with  commercial  airlines, 
the  military  services,  and  the  Civil  Aeronautics  Administration. 


METHODS  OF  GUST -STRUCTURE  MEASUREMENTS 


The  general  method  followed  in  the  investigation  of  gust  structure 
has  been  to  fly  airplanes  in  rough  air  and  to  deduce  the  gust  charac¬ 
teristics  from  the  reactions  of  the  airplane.  Whenever  possible,  the 
airplane  which  was  the  primary  instrument,  as  indicated  in  reference  1, 
has  been  calibrated  by  testing  dynamically  similar  scaled  models  in  the 
Langley  gust  tunnel  in  order  to  obtain  the  actual  relation  between 
acceleration  and  gust  intensity  for  various  gust  shapes.  Most  of  the 
data  on  over-all  gust  structure  have  been  obtained  from  records  of  accel¬ 
eration  and  airspeed  as  a  function  of  time-  These  records  have  provided 
information  on  the  gust  velocities,  the  gradient  distance,  and  the  gust 
spacing.  In  the  evaluation  of  the  records,  the  pitching  motion  of  the 
airplane  is  neglected  except  in  those  cases  where  it  can  be  included 
through  calibrations  made  in  the  gust  tunnel. 
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In  addition  to  measurements  of  over-all  gust  structure,  some 
measurements  have  teen  made  of  the  distribution  of  gust  velocity  along 
the  airplane  span  ty  recording  the  local  wing  pressures  at  various 
stations  along  the  span  and  the  airspeed.  From  thin-airfoil  theory,  the 
change  in  pressure  at  any  location  on  an  airfoil  divided  ty  the  dynamic 
pressure  q  is  a  function  of  the  angle-of-attack  change  and  that,  in 
turn,  is  equal  to  the  effective  gust  velocity  divided  ty  the  equivalent 
airspeed  of  the  airplane.  Local  values  of  true  gust  velocity  cannot  te 
measured  because  the  motions  of  the  airplane  are  somewhat  involved  and 
cannot  te  taken  into  account. 

The  intensity  of  horizontal  gusts  has  teen  obtained  from  airspeed 
records  on  the  assumption  that  the  absolute  velocity  of  the  airplane 
remains  constant  in  spite  of  rapid  changes  in  wind  speed  which  the  air¬ 
plane  may  experience.  The  apparent  changes  in  airspeed  are  taken  as  the 
horizontal  gust  velocities  when  not  associated  with  large  normal 
accelerations. 


APPARATUS  AND  TESTS 


Five  airplanes  having  the  characteristics  shown  in  table  I  and 
figure  6  have  teen  used  to  obtain  gust -structure  data.  As  may  te  noted 
from  figure  6,  all  airplanes  except  the  XBM-1  were  monoplanes,  and  the 
F-61C  airplane  had  twin  tail  booms  instead  of  the  usual  fuselage.  Table  I 
shows  that  the  weight  of  the  airplanes  varied  from  about  700  pounds 
to  55,000  pounds,  and  the  wing  mean  geometric  chord  varied  from  about 
4  feet  for  the  Aeronca  C-2  airplane  to  about  19  feet  for  the  XB-15  air¬ 
plane-  The  "sea-level"  mass  parameters  of  the  various  airplanes  varied 
from  about  6  to  23.  In  the  case  of  the  Aeronca  C-2  and  the  XBM-1  air¬ 
planes,  the  response  to  known  gusts  had  been  obtained  by  tests  in  the 
Langley  gust  tunnel  (references  7  and  8,  respectively) . 

For  the  determination  of  the  gust  intensity,  gradient  distance,  and 
gust  spacing,  all  airplanes  carried  a  standard  set  of  instruments 
consisting  of: 

(1)  Recording  accelerometer  located  at  the  center  of  gravity 

(2)  Airspeed  altitude  recorder 

(3)  Synchronizing  timer 
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The  instrument  ranges  were  adjusted  to  suit  the  particular  airplane. 

All  instruments  were  equipped  to  give  sufficient  record  time  and  were 
operated  at  film  speeds  to  give  adequate  time  resolution.  A  switch 
convenient  to  either  the  pilot  or  the  observer  enabled  the  instruments 
to  be  turned  on  or  off  at  will. 

All  airplanes  carried,  in  addition  to  the  instruments  mentioned, 
special  instruments  suited  to  the  particular  project.  The  203-15  air¬ 
plane,  for  example,  was  equipped  with  strain  gages  on  wing  beams,  and 
optigraphs  were  provided  to  measure  both  wing  and  tail  deflections. 

The  XC-35  airplane  carried  four  pressure  recording  units  that  were 
connected  to  orifices  located  on  the  wing  of  the  airplane  (see  fig.  7) 
to  measure  the  differential  pressure  between  the  upper  and  lower  sur¬ 
faces.  Radiosonde  data,  as  well  as  temperature  measurements,  were  also 
obtained  during  the  tests  of  the  XC~35>  XBM-1,  Aeronca  C-2,  and  F-61C 
airplanes.  Table  II  gives  the  scope  of  the  flight  tests-  With  the 
exception  of  the  flight  investigations  with  the  XB-15  airplane,  all  tests 
had  objectives  associated  with  research  on  gust  structure,  and  the  flight 
tests  were  under  the  control  of  the  test  engineer  as  to  the  performance 
and  the  type  of  flying  done.  The  stability  of  the  various  airplanes  was 
adjusted  whenever  possible  so  that  the  airplanes  showed  stick-fixed 
stability. 


RESULTS 


The  time  histories  of  acceleration  and  airspeed  obtained  during 
these  flight  investigations  were  evaluated  to  obtain  the  effective 
gust  velocity,  the  gradient  distance,  and  the  corresponding  true  velocity 
As  noted  previously,  every  acceleration  peak  can  be  evaluated  to  obtain 
the  effective  gust  velocity  corresponding  to  the  acceleration  increment 
as  measured  from  the  1  g  datum,  but  only  those  acceleration  peaks  that 
are  preceded  by  a  smooth  part  of  record  can  be  evaluated  to  obtain  the 
gust -gradient  distance  and  the  true  gust  velocity. 

Gust  intensity. -  The  available  data  on  gust  intensity,  as  measured 
by  the  effective  gust  velocities,  have  been  examined  in  considerable 
detail  in  reference  9,  and.  the  frequency  distributions  given  in  figure  8 
were  derived  therein.  The  results  given  in  figure  8  are  based  on  early 
German  work,  NACA  flight  tests,  and  V-G  records  from  commercial  transport 
operations.  Curves  A  and  B  are  the  approximate  limits  of  the  frequency 
distributions  of  effective  gust  velocities  for  over-all  operating  con¬ 
ditions  and  show  the  probability  that  a  gust  will  exceed  any  selected 
value  of  intensity. 

In  addition,  a  more  detailed  investigation  of  limited  scope  has 
been  made  (reference  10)  of  the  XC-35  data  to  compare  the  frequency 
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distributions  of  effective  gust  velocity  in  thunderstorms  and  line 
squalls  for  various  altitude  ranges.  Table  III  summarizes  the  data 
obtained  to  a  threshold  of  4  feet  per  second  for  several  ranges  of  alti¬ 
tude  from  the  surface  to  35*000  feet-  Table  III  also  lists  the  record 
time  and  the  average  spacing  between  gusts  for  each  altitude  range. 

Figure  9  shows  the  probability  distributions  of  effective  gust  velocity 
for  the  several  altitude  ranges  that  are  obtained  by  curves  fitted  to 
the  data  of  table  III.  In  this  investigation,  the  data  for  altitudes 
below  5000  feet  were  not  used  because  the  conditions  investigated  repre¬ 
sented  random  record  taking  in  both  clear  air  and  clouds  at  the  end  of 
the  flights. 

From  the  tests  performed  with  the  XC-35  and  the  F-61C  airplanes, 
the  variation  of  expected  maximum  effective  gust  velocity  with  altitude 
within  convective  clouds  has  been  determined.  These  results  are  sum¬ 
marized  in  figure  10,  in  which  the  data  for  the  XC-35  airplane  corre¬ 
spond  to  100  miles  of  rough-air  flying,  while  those  for  the  F-61C  air¬ 
plane  correspond  to  1000  miles.  The  latter  data  were  obtained  during 
the  194-6  operations  of  the  U.S.  Weather  Bureau  thunderstorm  project  at 
Orlando,  Fla.,  and  are  based  on  the  maximum  values  of  effective  gust 
velocity  for  each  3000-foot  interval;  however,  the  XC-35  data  utilized 
a  complete  count  of  all  gusts,. 

Gust  spacing. -  Some  information  on  gust  spacing,  that  is,  the 
distance  from  peak  to  peak,  has  also  been  obtained.  In  reference  11, 
an  analysis  was  made  of  the  XC-35  data  to  determine  the  average  spacing 
between  the  large  gusts  (defined  as  Ue  greater  than  5  to  8  fps) 
in  areas  of  continuous  rough  air.  Continuous  rough  air  was  defined 
as  a  sequence  of  large  gusts  in  which  the  spacing  between  gusts 
was  less  than  2.2  seconds.  Table  IV,  reproduced  from  reference  11, 
summarizes  some  of  the  pertinent  characteristics  of  the  large  isolated 
and  repeated  gusts.  The  table  presents  a  comparison  of  the  average  and 
the  range  of  intensities  and  the  average  and  range  of  spacings  for 
isolated  single  gusts  and  sets  of  two  and  three  repeated  gusts  occurring 
with  equal  frequency. 

More  general  information  on  gust  spacing,  contained  In  reference  9* 
indicates  that  the  average  spacing  between  gusts  that  were  counted  to  an 
estimated  threshold  of  0-3  foot  per  second  was  a  function  of  the  airplane 
mean  geometric  chord.  This  count  Indicated  an  average  spacing  between 
successive  peaks  of  11  chords  so  that  the  number  of  gusts  per  mile  of 
rough  air  is  roughly  equal  to  500/c.  The  definition  of  the  number  of 
gusts  per  mile  of  rough  air  led  to  the  concept  of  path  ratio  which 
provides  a  measure  of  the  percent  of  rough  air  encountered  under  oper¬ 
ating  conditions. 

Gust -gradient  distance. -  Figure  11  presents  a  typical  plot  of 
measurements  of  gradient  distance  and  gust  velocity.  Inspection  of  the 
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figure  indicates  that  the  scattered  pattern  of  data  is  not  amenable  to 
detailed  analysis  and,  for  such  purposes,  the  data  on  gust-gradient 
distance  were  sorted  according  to  the  variable  under  consideration. 
Figures  12(a),  12(b),  and  12(c)  present  the  data  on  the  average  gust- 
gradient  distance  as  a  function  of  the  gust  velocity  (width  of  bracket, 

4  fps) .  Figure  12(a)  is  a  summary  of  the  data  for  all  airplanes 
with  HaY  expressed  as  multiples  of  the  mean  geometric  chord} 

figure  12(b)  presents  the  same  data  with  HaY  expressed  in  multiples 
of  the  airplane  wing  span}  and  figure  12(c),  in  feet*  For  comparison 
with  other  figures  and  for  subsequent  use,  a  faired  curve  was  drawn 
through  the  data  in  figure  12(a)  to  agree  best  with  the  data  from 
the  XC- 35  airplane. 

The  results  presented  in  figures  12(a),  12(b),  and  12(c)  do  not 
include  the  complete  data  for  every  airplane  since  some  brackets,  parti¬ 
cularly  below  about  8  feet  per  second,  are  not  complete.  The  magnitude 
of  the  acceleration  increments  for  values  of  U  less  than  8  feet  per 
second  was  about  the'  same  order  as  the  resolution  of  the  accelerometer. 
For  the  larger  gust  velocities,  seme  "average"  values  were  not  plotted 
since  the  number  of  points  represented  was  less  than  three. 

In  some  cases,  sufficient  data  were  available  to  obtain  the  probable 
values  of  the  gust-gradient  distance.  The  probable  value  was  determined 
by  fitting  a  theoretical  frequency  distribution  to  the  data  and  deter¬ 
mining  the  point  with  the  highest  frequency.  The  results  are  plotted 
against  the  corresponding  average  gust -gradient  distance  HaT  in 
figure  13.  Figure  13  includes  two  solid  lines,  a  45°  line  representing 
the  equality  of  the  quantities  and  a  line  through  the  data  offset  by 
2  chords . 

Data  for  the  gust-gradient  distance  for  the  XC-35  and  F-61C  air¬ 
planes  have  been  utilized  to  determine  the  variations  of  the  average 
gradient  distance  expressed  in  chords  as  a  function  of  altitude.  (See 
fig.  14.)  In  the  case  of  the  XC-35  airplane  where  no  fixed  altitudes 
were  flown,  altitude  brackets  were  selected,  and  the  average  gradient 
distance  was  plotted  against  the  center  value  for  the  pertinent  bracket. 
The  F-61C  airplane  tests  were  made  at  fixed  altitudes,  and  the  data  have 
been  plotted  at  the  appropriate  altitude  level. 

Spanwise  guet  distribution. -  Records  of  local  wing  pressures  and 
airspeeds  obtained  during  the  flight  investigation  with  the  XC-35  air¬ 
plane  (reference  12)  were  evaluated  to  obtain  the  local  values  of 
effective  gust  velocity  at  each  of  four  spanwise  stations.  The  general 
method  followed  was  to  plot  the  data  for  each  gust  as  a  function  of  the 
spanwise  location  and  connect  the  points  by  straight  lines.  The  indivi¬ 
dual  spanwise  distributions  were  later  identified  with  the  corresponding 
values  of  the  gust  velocity  U  and  the  effective  gust  velocity  Ue 
obtained  from  the  other  flight  records. 
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Data  on  the  spanwise  distribution  were  sorted  according  to  the  shape 
of  the  distribution  as  indicated  in  figure  15*  Six  shapes  ranging  from 
the  uniform  or  rectangular  distribution  to  the  double  triangular  distri¬ 
bution  are  shown,  and  under  each  shape  is  noted  the  total  number  of  gusts 
classified  in  that  category.  Figure  15  also  indicates  possible  variations 
included  under  each  heading. 

The  spanwise  gust  distributions  represented  in  figure  15  were,  in 
turn,  evaluated  graphically  to  obtain  values  of  the  lateral  gust -gradient 
distance  H.  The  results  are  shown  in  figure  1 6  where  the  average  values 
of  the  gradient  distance  in  chords  are  given  for  a  range  of  values  of 
gust  velocity.  The  associated  gust  velocity  U  was  obtained  from  the 
synchronized  acceleration  data,  and  the  faired  curve  shown  in  figure  1 6 
corresponds  to  that  given  in  figure  12(a) . 

The  significance  of  the  values  of  Ue  obtained  from  pressure 
measurements  as  compared  to  those  determined  from  acceleration  measure¬ 
ments  has  been  examined  in  reference  12.  Figure  17,  reproduced  from 
that  reference,  shows  the  average  value  of  the  local  effective  gust 
velocity  as  a  function  of  the  effective  gust  velocity  for  a  given  gust 
as  determined  from  acceleration  records-  The  average  value  of  Ue  from 
the  pressure  measurements  is  the  average  value  weighted  according  to  the 
amount  of  wing  area  represented  by  the  particular  orifice  or  measuring 
station.  Also  -shown  in  the  figure  is  the  line  for  equality  of  the  two 
measures  of  Id  and  the  error  band. 

Since,  in  some  cases,  it  is  convenient  to  divide  the  lateral  gust 
distribution  into  symmetrical  and  unsymmetrical  components,  the  gust 
shapes  listed  in  figure  15  were  classified  according  to  whether  they  were 
symmetrical  or  unsymmetrical  and  then  represented  as  having  a  linear 
variation  of  gust  velocity  across  the  span  of  the  airplane-  The  results 
are  shown  in  figure  18.  The  inset  in  figure  18  indicates  the  type  of 
variation  assumed  and  the  two  quantities  obtained  from  the  evaluation. 

The  average  gust  velocity  U„  and  the  increment  in  gust  velocity  at 

°av 

each  wing  tip  AUe  were  taken  as  the  symmetrical  and  unsymmetrical 
components,  respectively.  Many  small  values  have  been  omitted  from 
the  figure. 

As  in  the  case  of  the  gust -gradient  distance,  the  data  have  been 
further  analyzed  to  find  an  average  unsymmetrical  gust  increment. 

Figure  19  has  been  prepared  from  the  summary  of  data  for  the  different 
gust  shapes  to  show  the  average  gust  velocity  for  a  gust  of  the  unsym¬ 
metrical  shape  as  a  function  of  the  average  gust  velocity  of  a  triangular 
gust  for  equal  frequencies  of  occurrence. 

In  addition  to  the  determination  of  spanwise  gust  distribution  from 
local  pressure  distributions,  synchronized  records  from  accelerometers 
located  at  the  center  of  gravity  of  the  airplane  and  in  the  wing  were 
.evaluated  to  obtain  the  angular  acceleration  in  roll,  which  is  related 
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to  the  spanwise  gust  distribution.  The  results  of  the  evaluation  of 
accelerometer  records  for  17  flights  are  shown  in  table  V,  which  indi¬ 
cates  the  frequency  of  occurrence  of  given  values  of  angular  accelera¬ 
tion  with  the  associated  values  of  normal  acceleration  increment. 

Longitudinal  gusts. -  Rapid  airspeed  fluctuations  have  been  used  to 
evaluate  the  magnitude  of  longitudinal  gusts.  Figure '20,  prepared  from 
an  analysis  of  the  XC-35  airplane  data,  compares  the  maximum  values  of 
the  horizontal  gust  velocities  obtained  from  the  rapid  airspeed  fluctua¬ 
tions  with  the  maximum  vertical  gust  velocities.  Each  point  represents 
a  separate  traverse  or  rim.  The  results  of  a  detailed  analysis  of  air¬ 
speed  and  acceleration  records  in  gusty  air  from  the  XBM-1  and  Aeronca  C-2 
airplane  investigations  are  shown  in  figure  21  and  indicate,  for  both 
airplanes,  the  relation  of  the  horizontal  gust  velocities  to  the  vertical 
gust  velocities  for  equal  frequencies  of  occurrence. 


ACCURACY  OF  RESULTS 


Consideration  of  instrumental  and  reading  errors  in  evaluating  the 
data  from  the  several  airplanes,  together  with  a  knowledge  of  the  problems 
involved  in  the  reactions  of  an  airplane,  indicates  that  the  possible 
errors  in  the  derived  quantities  that  are  of  interest  are  approximately 
as  follows: 


Gust  velocity  based  on  acceleration,  percent  . .  ±10 

Gust  velocity  from  pressure  measurements,  feet 

per  second .  ±4 

Gust  velocity  from  airspeed  measurements,  feet 

per  second . . . ±10 

Gust -gradient  distance,  chords . ±1  to  ±2 

Lateral  gust-gradient  distance,  percent  .  ±20 


Spacing,  chords 


The  accuracy  of  the  gust -velocity  measurements,  whether  true  gust 
velocity  or  effective  gust  velocity,  from  acceleration  records  depends 
on  the  quality  of  the  records  obtained,  the  availability  of  a  gust- 
tunnel  calibration  of  the  airplane,  and  experience  in  reading  and 
interpreting  the  records. 


The  local  gust  velocities  obtained  from  wing-pressure  measurements 
were  tested  for  accuracy  by  comparing  the  average  value  of  the  local 
gust  velocities  with  the  effective  gust  velocity  from  acceleration  data. 
The  results  indicate  that  the  two  quantities  are  equivalent  but  that  the 
data  scatter  over  a  range  of  about  4  feet  per  second. 
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The  gradient  distances  determined  from  accelerometer  records  are  a 
function  of  loth  the  time  resolution  and  the  response  characteristics  of 
an  airplane  due  to  known  gusts.  Examination  of  gust-tunnel  test  results 
indicates  that,  under  average  conditions,  the  airplane  response  places 
two  limits  on  the  gust-gradient  distance  determined  fhom  acceleration 
records.  The  first  limit,  for  the  shorter  gradient  distances,  is  speci¬ 
fied  by  the  lag  in  lift.  The  gradient  distances  lying  above  the  dash 
line  in  figure  4  are  not  generally  recognizable  on  an  accelerometer 
record.  Thus,  this  limit  is  from  2  chords  to  about  4  chords  for  the 
average  airplane.  The  second  limit,  for  the  longer  gradient  distances, 
arises  from  the  fact  that,  as  the  airplane  penetrates  farther  and  farther 
into  a  long  gradient  gust,  the  effect  of  the  pitching  and  vertical  motion 
of  the  airplane,  under  the  action  of  the  gust,  tends  to  counteract  the 
contribution  of  lift  due  to  each  increment  of  gust  as  it  is  encountered. 

The  lateral  gust -gradient  distance  is  a  more  questionable  quantity 
than  the  gust-gradient  distance  determined  from  accelerometer  records 
since  no  single  concept  as  to  shape  has  been  obtained.  The  lack  of 
detailed  data  and  the  erratic  character  of  the  gust  distributions  in  a 
spanwise  direction  (fig-  15)  require  the  use  of  considerable  Judgment 
in  arriving  at  any  numerical  values.  Comparison  of  the  actual  values  of 
the  gust-gradient  distance  from  accelerometer  records  and  the  corre¬ 
sponding  values  of  the  lateral  distance  from  pressure  records  indicates 
that  the  data  are  reasonable  and  in  agreement  with  the  actual  conditions. 
From  the  preceding  remarks  it  is  obvious  that  the  error  is  significant 
and  the  value  of  20  percent  previously  mentioned  is  only  a  crude  estimate* 


DISCUSSION 


Gust  intensity. -  Because  only  a  small  number  of  all  the  acceleration 
peaks  in  rough  air  can  be  evaluated  for  true  gust  velocities,  the  effective 
gust  velocity  has  been  used  as  the  measure  of  intensity  and  frequency  of 
atmospheric  gusts.  Although  the  effective  gust  velocity  Ue  is  a 
fictitious  quantity,  it  bears  a  fixed  relation  to  the  gust  velocity  for 
a  given  airplane  and  given  conditions  of  air  density,  gust  shape,  and 
gust  size.  The  effective  gust  velocity  is  generally  about  50  percont  to 
70  percent  of  the  gust  velocity. 

The  largest  effective  gust  velocity  recorded  to  date  was  about 
55  feet  per  second  and  corresponds  to  a  gust  velocity  of  100  feet  per 
second  for  typical  conditions.  Fortunately,  gust  velocities  this  large 
are  not  frequently  encountered.  The  distributions  of  effective  gust 
velocity  in  figure  8  show  that  the  limiting  distributions  for  a  wide 
range  of  airplane  sizes  and  weather  conditions  are  reasonably  close. 
Deference  9  concluded,  on  the  basis  of  these  results,  that  the  distri¬ 
bution  of  gust  intensities  was  essentially  independent  of  airplane  size 
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and  source  of  turbulence.  Figure  9,  in  contrast,  shows  that  the  XC-35 
data  for  different  altitude  brackets  result  in  distributions  with  slopes 
that  vary  and,  when  checked  against  figure  8,  differ  from  those  curves. 

This  discrepancy  indicates  that  the  distributions  of  effective  gust 
velocity  may  not  be  independent  of  the  source  of  turbulence.  Independ¬ 
ency  is  a  reasonable  assumption  at  this  time  for  general  use. 

Figures  9  and  10  show  that,  for  altitudes  above  9000  feet,  the 
maximum  gust  intensity  in  rough  air  associated  with  convective  clouds  is 
essentially  a  constant.  Below  9000  feet,  the  XC-35  data  indicate  a 
decrease  in  intensity,  due  in  part,  according  to  reference  10,  to  the 
fact  that  the  records  were  taken  in  clear  air  and  clouds.  Consideration 
of  figure  10  indicates  that  little  difference  exists  in  the  expected 
maximum  gust  velocities  for  all  altitude  ranges. 

Gust -gradient  distance . -  Figure  11  and  similar  plots  for  other 
airplanes  show  that  gust-gradient-distance  data  are  of  a  random  character. 
Comparisons  of  such  plots  indicate  that  airplane  size  might  be  a  signifi¬ 
cant  parameter,  as  in  the  case  of  boats  where  waves  of  small  wave  length 
that  are  of  no  concern  to  a  large  boat  cause  the  small  boat  considerable 
difficulty.  In  a  similar  manner,  gusts  that  cause  a  rough  ride  on  a 
small  airplane,  such  as  the  Aeronca  C-2,  would  be  expected  to  be  of  such 
small  size  as  to  have  little  effect  on  a  large  airplane,  such  as  the  XB-15- 

Comparison  of  figures  12(a),  12(b),  and  12(c)  indicates  that  the 
gradient-distance  data  show  much  less  scatter  on  the  basis  of  the  mean 
geometric  chord  than  on  the  basis  of  span  or  feet.  The  scatter  of  points 

when  the  gust-gradient  distance  is  defined  in  chords  is  about  l^  to  1 

for  a  given  value  of  U,  while  for  the  gust-gradient  distance  plotted  in 
terms  of  airplane  span,  the  scatter  is  greater.  The  data  of  figure  12(c) 
show  that  when  the  gust-gradient  distance  is  expressed  in  feet,  the 
average  values  of  the  distance  vary  from  40  feet  to  200  feet  for  a  given 
gust  intensity.  Inasmuch  as  the  agreement  for  the  different  sets  of 
data  is  best  when  the  gradient  distance  is  expressed  in  chords,  the  gust- 
gradient  distance  appears  to  be  roughly  independent  of  the  airplane  when 
expressed  in  wing  mean  geometric  chords. 

The  adequacy  of  the  average  value  of  the  gust-gradient  distance  for 
the  representation  of  gradient -distance  data  deserves  some  consideration. 
The  relation  between  the  average  value  and  the  most  probable  value  shown 
in  figure  13  indicates  that  the  probable  gradient  distance  is  about  1  chord 
to  2  chords  less  than  the  average  gradient  distance.  Since  the  offset 
appears  to  be  constant,  the  use  of  the  average  gust-gradient  distance, 
which  is  generally  easier  to  determine,  is  believed  to  be  satisfactory 
for  analysis. 
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As  Indicated  in  tail e  II,  a  variety  of  weather  conditions,  from  line 
squalls  to  gustiness  in  the  ground  houndary  layer,  is  included  in  the 
present  data.  If  the  data  of  figure  12(a)  are  assumed  to  he  comparahle, 
regardless  of  airplane,  the  type  of  weather  appears  to  have  no  effect  on 
the  gust -gradient  distance. 

Figure  14  indicates  that  the  average  gust-gradient  distance  is  inde¬ 
pendent  of  altitude  in  convective  clouds.  The  data  scatter  somewhat  hut 
the  over-all  agreement  appears  to  he  good. 

General  considerations  would  indicate  that  the  gradient  distance 
might  increase  with  gust  velocity.  Previous  study  and  figure  12(a)  have 
indicated  such  a  relation,  although  the  evidence  is  not  conclusive. 

When  the  faired  curves  of  figures  12(a)  and  16  are  considered,  such  a 
variation  appears  to  he  a  reasonable  estimate-  The  evidence  available 
(fig.  12(a))  is  believed  to  indicate  that  for  large  gusts  the  average 
gradient  distance  is  essentially  constant  hut  decreases  rapidly  at  the 
smaller  gust  velocities. 

S uanwise  gust  distribution. -  The  shapes  of  the  lateral  gust  distri¬ 
butions,  as  indicated  by  figure  15,  take  many  irregular  forms.  From 
consideration  of  the  frequency  of  occurrence,  the  triangular,  double- 
triangle,  and  unsymmetrical  gust  shapes  predominate.  Inasmuch  as  the 
double  triangle  can  lead  to  wing  bending  moments  along  the  span  equal 
to  those  due  to  a  uniform  gust  distribution,  while  the  triangular  gust 
shape  would  lead  to  reduced  wing  bending  moments,  the  selection  for 
design  of  a  uniform  gust  velocity  across  the  span  for  the  symmetrical 
load  condition  is  conservative.  The  relative  percentages  of  the  sym¬ 
metrical  (triangle,  and  so  forth)  and  unsymmetrical  gust  shapes  indicate 
that  the  unsymmetrical  gust  varying  uniformly  across  the  span  appears 
quite  frequently. 

Comparison  of  the  data  for  the  average  lateral  gust -gradient  dis¬ 
tance  as  a  function  of  gust  velocity  U  (fig-  16)  with  corresponding 
results  for  longitudinal  gust-gradient  distances  (fig-  12(a))  Indicates 
that  the  gradient  distances  are  essentially  the  same  in  both  directions. 
This  agreement  is  indicated  by  the  fact  that  the  same  empirical  curve 
fits  both  sets  of  data  and  the  results  show  the  same  trend  of  increasing 
gust -gradient  distance  with  increasing  gust  velocity.  Because  data  from 
only  one  airplane  have  been  presented  and  the  precision  of  measurement 
of  lateral  values  Is  poor,  aB  previously  cited,  the  excellent  agreement 
between  figures  16  and  12(a)  is  to  same  degree  fictitious,  the  degree 
being  unknown.  Within  the  limitations  of  results  the  lateral  and  longi¬ 
tudinal  gust  dimensions  are  concluded  to  be  essentially  the  same- 

Inspection  of  figure  18,  relating  the  unsymmetrical  components  of 
the  effective  gust  velocity  to  the  average  uniform  gust  velocity  across 
the  span,  indicates  that  no  clear  relations  between  the  unsymmetrical 
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components  and  the  average  gust  velocity  exists  "but  that  two  "boundaries 
appear  to  exist.  A  fairly  definite  "boundary  rising  from  the  lower  left- 
hand  corner  of  the  figure  and  corresponding  to  equality  of  the  unsym- 
metrical  component  of  gust  velocity  and  the  average  gust  velocity  is 
indicated.  As  higher  values  of  the  average  gust  velocity  are  considered, 
the  results  scatter  considerably,  but  an  upper  boundary  is  found  at 
from  9  to  10  feet  per  second.  The  boundary  holds  for  gust  velocities 
up  to  40  to  48  feet  per  second.  On  the  basis  of  this  upper  limit,  the 
unsymmetrieal  component  of  gust  velocity  above  10  feet  per  second  can  be 
considered  rare-  Therefore,  the  unsymmetrieal  gust  can  be  considered  to 
be  composed  of  two  components,  an  unsymmetrieal  component  amounting  to 
a  maximum  of  about  >10  feet  per  second  with  opposite  signs  imposed  at  each 
wing  tip  and  a  linear  gradient  across  the  span,  and  an  average  uniform 
gust  velocity  of  any  pertinent  value. 

A  question  of  some  importance  is  the  determination  of  the  magnitude 
of  the  unsymmetrieal  components  for  other  aircraft  of  entirely  different 
spans.  For  example,  if  the  span  of  the  airplane  were  doubled,  would  the 
unsymmetrieal  component  amount  to  twice  the  value  obtained  from  the  data 
shown  herein?  In  order  to  obtain  information  on  this  question,  the 
maximum  value  of  the  unsymmetrieal  gust  components  computed  in  refer¬ 
ence  13  from  angular -acceleration  data  was  compared  with  the  values 
obtained  herein.  It  was  found  that  for  the  XB-15  airplane,  which  has  a 
span  of  about  150  feet,  the  unsymmetrieal  components  of  gust  velocity 
at  the  wing  tip  would  be  about  13-7  feet  per  second.  This  value  is 
slightly  greater  than  indicated  by  figure  l8  but  is  not,  by  any  means, 
proportional  to  the  span  of  the  airplane.  Similar  data  of  a  statistical 
nature  obtained  for  the  XC-35  airplane  (table  V)  indicate  that  the  maxi¬ 
mum  value  of  the  angular  acceleration  on  that  airplane  varied  from  about 
1.95  "to  about  2.25  radians  per  second  per  second.  Computations  utilizing 
the  formula  given  in  reference  13  yield  a  maximum  gust -velocity  component 

of  about  12?j  feet  per  second.  This  value  is  in  fair  agreement  with  the 

one  obtained  for  the  XB-15  airplane  and  indicates  that  the  tip  values 
are  independent  of  span.  Since  the  values  based  on  angular  accelera¬ 
tion  are  somewhat  more  conservative  than  those  based  on  wing  pressures, 
more  reliance  should  be  placed  on  these  values. 

"When  the  information  on  unsymmetrieal  gusts  is  utilized,  the  average 
effective  gust  velocity  that  should  be  used  in  conjunction  with  the  maxi¬ 
mum  unsymmetrieal  component  of  14  feet  per  second  must  be  determined. 

The  data  of  figure  19  indicate  that  for  an  effective  gust  velocity  in  the 
neighborhood  of  30  feet  per  second,  the  average  gust  velocity  for  the 
unsymmetrieal  gust  would  be  about  22  feet  per  second. 

Gust  spacing. -  Inspection  of  table  IV  indicates  that  the  spacing  of 
repeated  gusts  for  two  or  three  gusts  in  sequence  is  about  22  chords 
with  a  spread  in  actual  values  of  some  20  chords  for  all  sequences. 
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The  spacing  is  approximately  twice  the  average  gust -gradient  distance  of 
10  to  l4  chords.  These  data  indicate  that  the  gust  shape  in  the  direc¬ 
tion  of  flight  is  either  triangular  or  sinusoidal  in  character  as 
contrasted  to  the  flat -top"  gust  assumed  in  past  years. 

From  time  to  time,  comparisons  are  made  of  the  gust  spacing  in 
reference  11  of  22  chords  with  the  gust  spacing  that  was  derived  from 
statistical  gust  data  in  reference  9  of  11  chords.  The  discrepancy  be¬ 
tween  the  two  figures  arises  from  the  fact  that  the  data  on  gust  spacing 
presented  in  reference  11  were  obtained  for  the  larger  gusts  and  represent 
data  on  gust  intensities  ranging  from  5  feet  per  second  up  to  the  maxi¬ 
mum  value  recorded,  whereas  the  gust  spacing  listed  in  reference  9 
is  based  on  the  average  spacing  of  all  gusts  from  a  gust  velocity  of 
0.  3  foot  per  second  to  the  maximum  value  experienced  by  any  of  the  air¬ 
planes  considered. 

The  sequences  of  gusts  found  in  reference  11  were  such  that  in 
about  two-thirds  of  the  cases  of  sets  of  two  the  gusts  were  of  opposite 
sign  while  the  remaining  sets  of  two  were  composed  of  gusts  of  like  sign. 
The  larger  gust  intensity  could  be  found  at  any  point  in  the  sequence 
under  consideration. 

Longitudinal  gusts.  -  Inspection  of  the  data  for  the  maximum  values 
of  the  horizontal  and  vertical  gust  velocities  in  the  same  traverse 
(fig.  20)  indicates  equality  of  the  gust  intensities  in  the  two  direc¬ 
tions.  Figure  21  for  the  XBM-1  and  Aeronca  C-2  airplanes  gives  essen¬ 
tially  the  same  indication,  but  the  agreement  is  not  exact  and  the 
discrepancy  between  the  values  of  the  gust  velocities  is  an  offset  of 
roughly  1  foot  per  second.  The  data  obtained  from  the  Aeronca  C-2  air¬ 
plane  show  the  larger  discrepancy.  Some  of  this  discrepancy  may  be  due 
to  the  fact  that  the  amount  of  horizontal-gust  data  was  larger  than  the 
amount  of  vertical -gust  data  for  the  Aeronca  C-2  airplane  and  smaller 
than  the  amount  for  the  XBM-1  airplane. 

Since  the  maximum  gust  intensity  in  the  horizontal  direction  is 
equal  to  the  vertical  gust  intensity  for  the  same  region  of  turbulence 
and  a  preliminary  investigation  of  frequency  distributions  indicates 
that  the  frequency  distributions  are  essentially  the  same,  the  atmos¬ 
pheric  turbulence  seems  to  be  isotropic.  Previous  indications  as  to 
structure  of  vertical  gusts,  that  is,  gust  spacing,  gust-gradient 
distance,  and  lateral  distributions  of  gust  velocity,  would  be  assumed 
to  apply  to  horizontal  gusts. 

The  preceding  results  indicate  equal  gradient  distances  for  the 
longitudinal  and  lateral  faces  of  the  gust  velocity  distribution,  and 
this  equality  suggests  that  the  gust  is  symmetrical  about  a  vertical 
axis.  In  keeping  with  this  concept,  the  findings  suggest  that  the 
distribution  might  be  visualized  as  a  four- sided  pyramid  with  a  base 
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length  of  20  to  28  chords.  While  such  a  velocity  distribution  may 
represent  the  average  gust,  because  of  the  wide  variations  in  shape,  a 
wedge-shaped  velocity  distribution  such  that  the  velocity  distribution 
along  the  span  is  uniform  and  the  longitudinal  shape  is  either  triangular 
or  sinusoidal  is  recommended  for  load  calculations. 


CONCLUDING  REMARKS 


Although  much  of  the  information  on  gust  structure  can  be  ration¬ 
alized  to  obtain  a  "standard"  gust,  this  material  should  be  used  with 
caution  lor  unconventional  types  of  aircraft.  The  gust  structure  of 
interest  is,  of  course,  dependent  to  some  degree  on  the  airplane  charac¬ 
teristics.  Many  other  sizes  of  gusts  exist  and  those  that  affect  the 
airplane  may  be  only  a  small  part  of  the  random  variations  that  exist  in 
the  atmosphere • 

Available  information  on  the  structure  of  atmospheric  gusts  has  \ 

shown  that,  when  the  gust  size  and  the  gust-gradient  distance  are  expressed  \ 
in  mean  geometric  chords,  the  gust  size  is  independent  of  airplane,  -  j 

weather,  topography,  and  altitude.  The  probable  size  of  the  gust  is 
20  chords,  and  the  probable  gust-gradient  distance  for  the  standard  gust  ■ 
is  about  10  chords.  A  wedge-shaped  gust  with  the  gust  velocity  uniform 
across  the  span  and  either  triangular  or  sinusoidal  in  shape  with  a  base  J 
of  20  chords  is  believed  to  be  the  proper  type  for  most  load  calculations^ 

For  a  sequence  of  gusts;  the  gusts  may  be  of  either  like  or  unlike 
sign  and  will  be  continuous,  with  the  gust  peaks  spaced  20  chords  apart. 

Since  the  maximum  horizontal  and  vertical  gust  intensities  and 
frequency  distributions  are  essentially  the  same  within  any  region  of 
rough  air,  the  gust  structure  obtained  for  vertical  gusts  should  apply 
equally  well  for  horizontal  gusts.  Such  information  is  pertinent  where 
gust  loads  are  under  consideration  for  diving  airplanes  or  missiles. 

Since  the  data  are  influenced  by  the  reactions  of  the  airplane,  the 
application  to  unconventional  configurations  should  be  investigated  by 
a  detailed  analysis  of  several  possible  combinations  of  gust  shape  and 
size. 


AIRPLANE  REACTIONS 


The  second  phase  of  the  gust-load  problem  is  that  of  determining 
the  reaction  or  forces  imposed  on  an  airplane  due  to  a  known  gust.  The 
factors  considered  are  the  aerodynamic  coefficients  and  the  possible 
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effects  of  stability  and  elasticity.  The  purpose  of  this  section  is  to 
indicate  what  is  known  about  these  factors  by  coordinating  the  available 
information.  In  most  cases,  specific  procedures  or  numbers  to  be  used 
in  load  calculations  are  not  obtained.  In  regard  to  some  factors,  the 
state  of  available  information  is  unsatisfactory  but  the  factor  is 
considered  and  its  status  is  indicated. 


METHODS 


The  three  methods  utilized  in  the  study  of  airplane  reactions  were  an 
analysis  and  two  experimental  methods.  The  experimental  methods  consist 
of  tests  of  models  in  the  Langley  gust  tunnel  and  flight  tests  with  full- 
scale  airplanes. 


Analysis 

Inspection  of  the  equations  given  in  references  lk  and  15  indicates 
that  the  inclusion  of  unsteady  lift  leads  to  variable  coefficients,  which 
result  in  integral  equations  when  the  unknown  appears  under  the  integral 
sign.  In  all  cases  the  integrals  that  appear  are  of  the  same  form  and 
represent  the  application  of  the  principle  of  superposition  to  unit- Jump 
solutions  to  obtain  the  response  to  arbitrary  or  known  disturbances. 

The  principle  cited  is  illustrated  in  figure  22  for  a  linear  varia¬ 
tion  in  angle  of  attack  for  which  the  lift  on  an  airfoil  after  s-j_  chords 
of  penetration  is  desired.  In  figure  22(a)  the  angle  of  attack  is  assumed 
to  vary  directly  with  s  and  the  development  of  lift  per  unit  change  in 
angle  of  attack  is  assumed  similar  to  that  shown  in  figure  3.  The  angle- 
of-attack  variation  is  assumed  to  be  approximated  by  a  series  of  unit 
changes  in  angle  of  attack  superimposed  in  the  s  direction.  For  steady 
lift,  the  corresponding  variations  of  lift  with  s  are  indicated  by  the 
straight  line  and  steps  in  figure  22(b).  The  unsteady  lift  develops  for 
each  step  in  angle  of  attack  according  to  the  dash  curves  and  the  approxi¬ 
mate  lift  at  sj_  is  the  sum  of  the  contributions  of  each  step  at  s^- 
The  result  at  s^  can  be  written  as 
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or  for  an  analytic  expression  of  differential  increments 


AC 
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da  , 
—  ds 
ds 


This  expression  is  commonly  known  as  Duhamel's  integral  and  is  illus¬ 
trated  by  Berg  in  reference  1 6.  The  integral  can  "be  evaluated  analy¬ 
tically  step  by  step  as  indicated  in  figure  22  or.  graphically  by  Carson's 
theorem  (reference  17) • 


As  previously  mentioned,  the  type  of  equation  obtained  is  of  the 

form 


s  r  ds 


ds 
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is  the  local  value  of  at  s-^.  The  difficulty  in 

solving  the  equation  arises  from  the  fact  that  the  relation  between  CL 

w 

and  s  must  be  known  before  the  second  integral  on  the  right  side  can 
be  evaluated.  The  equation  shown  is  a  simple  case  and  the  more  complete 
equations  of  motion  consist  of  additional  integrals  and,  in  some  cases, 
derivatives  of  integrals. 


where  Ct 

"1 


Solutions  of  equations  of  the  type  shown  have  been  obtained  by 
methods  of  successive  approximation,  Fredholm's  solution  (reference  6) , 
Laplace  transforms  (reference  14) ,  or  by  assuming,  as  in  reference  18, 
that  the  function  is  known.  The  need  for  solutions  for  arbitrary  dis¬ 
turbances  and  the  complicated  nature  of  solutions  generally  lead  to 
graphical  or  approximate  methods  in  all  but  a  few  cases. 


Gust -Tunnel  Testing 

The  Langley  gust  tunnel  (fig.  23)  was  built  to  permit  the  determina¬ 
tion  of  airplane  reactions  and  other  pertinent  quantities  under  controlled 
and  known  conditions.  It  consists  of  a  catapult  to  launch  a  dynamically 
scaled  airplane  model  into  steady  level  flight  through  a  vertical  jet  of 
air  having  characteristics  that  are  under  control,  a  means  of  catching 
the  model,  and,  finally,  suitable  equipment  to  record  the  required 
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information.  The  old  KACA  gust  tunnel  described  in  reference  19  was 
capable  of  handling  3“foot-span  models  of  airplanes  at  speeds  up  to 
about  50  miles  per  hour  but  in  19^5  was  replaced  by  a  tunnel  that  is 
able  to  handle  6-foot-span  models  at  speeds  up  to  100  miles  per  hour. 

The  accuracy  of  measurement  at  present  is  about  0.05g  for  accelera¬ 
tion,  0-3  foot  per  second  for  speed,  and  0.01  inch  for  deflection  of  the 
wing  or  flap.  Angular  motions  are  determined  within  0.1°  of  pitch-angle 
increment  and  about  0.2°  for  the  position  of  any  flap.  All  measurements 
except  that  of  speed  are  incremental  values  from  steady  conditions  and, 
as  such,  depend  upon  the  steadiness  of  flight  prior  to  entry  into  the 
gust.  At  peak  acceleration  for  sharp-edge  gusts,  the  acceleration  incre¬ 
ment  is  considered  accurate  to  within  O.lg  and  the  pitch-angle  increment 
to  about  0.1°,  but  the  accuracy  is  questionable  for  the  longest  gust- 
gradient  distance  of  about  16  chords. 

In  fundamental  studies  such  as  the  determination  of  unsteady-lift 
functions,  the  solution  obtained  from  gust-tunnel  tests  is  an  indirect 
one  since  the  acceleration  increment  due  to  a  gust  is  dependent  on  the 
difference  between  two  terms.  Thus,  it  is  not  possible  to  determine 
whether  the  correct  values  of  the  unsteady-lift  functions  are  used,  but 
only  whether  the  use  of  unsteady-lift  functions  for  a  given  shape  is 
adequate.  If  sufficient  data  are  collected  utilizing  different  models 
of  many  sizes  and  shapes  in  different  gusts,  the  various  unsteady -lift 
functions  can  be  evaluated  to  a  limited  degree. 


Flight  Investigations 

Two  approaches  have  been  utilized  for  flight  investigations  of  air¬ 
plane  reactions,  one  consisting  simply  of  detailed  analysis  of  flight 
records  to  obtain  the  relations  d.esired,  for  example,  between  the  pitch 
of  an  airplane  and  the  acceleration  imposed  on  the  airplane,  and  the 
other  consisting  of  flying  an  airplane  in  rough  air  to  obtain  statistical 
data  for  different  conditions  and  then  statistically  comparing  the  reac¬ 
tions  of  the  airplane.  The  first  approach  was  used  in  a  rather  crude 
manner  in  reference  3>  in  which  a  determination  was  made  of  the  effect 
of  airplane  wing  loading  and  speed  on  the  gust  load  factor.  The  statis¬ 
tical  approach  may  permit  checks  to  be  made  of  both  theoretical  and  gust- 
tunnel  results,  but  the  procedures  and  techniques  have  not  been  developed 
to  the  point  where  precise  results  can  be  obtained. 


TRANSIENT  AER OHfNAMICS 


The  material  presented  in  this  section  covers  the  development  of  the 
lift  force  and  downwash  but  does  not  cover  the  variation  in  moment 
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coefficient.  The  development  of  lift  following  a  sudden  change  in  angle 
of  attach  is  assumed  to  he  independent  of  the  slope  of  the  lift  curve 
and  can  he  treated  as  a  separate  subject.  The  neglect  of  possible  varia¬ 
tions  in  moment  coefficient  appears  justified  at  this  time  on  the  basis 
that  theory  indicates  that  the  zero-lift  moment  coefficient  is  unaffected 


Unsteady-Lift  Functions 

In  the  present  paper,  the  variation  of  lift  coefficient  for  a  unit 
change  in  angle  of  attack  has  been  referred  to  as  Ch  or  as  the 

ija 

Wagner  function,  and  the  corresponding  variation  in  lift  coefficient 
following  a  change  in  angle  of  attack  due  to  penetrating  a  unit  gust  has 

been  called  Ct  or  the  Kus sner  function.  The  Kiissner  function  is 
g 

related  to  the  Wagner  function  in  that  it  has  been  derived  from  it  by 
considering  the  airfoil  penetrating  a  gust  to  be  replaced  by  a  deforming 
airfoil  where  the  deformations  correspond  to  the  chordwise  angle -of - 
attack  distribution  due  to  gust  penetration.  Unless  otherwise  stated, 
the  lift  at  any  time  after  the  start  of  the  angle-of-attack  change  is 
expressed  as  a  fraction  of  the  final  lift  coefficient. 

Infinite  aspect  ratio.  -  Both  the  Wagner  function  and  the  Kiissner 
function  for  unsteady  lift  have  been  derived  a  number  of  times.  The 
former  was  originally  given  by  Wagner  in  reference  20  and  has  been 
checked  experimentally  by  Walker  in  reference  21  by  measuring  the  circula 
tion  about  a  wing  following  a  sudden  change  in  forward  speed-  Figure  24 
summarizes  the  computed  and  observed  circulations  obtained  by  Walker. 
Although  the  amount  of  direct  experimental  verification  is  limited, 
several  analytical  studies  agree  and  Indicate  that  the  Wagner  function 
should  be  close  to  correct. 

The  Kiissner  function  has  also  been  computed  a  number  of  times  and 
different  results  have  been  obtained.  Figure  25,  which  is  based  on  data 
taken  from  references  2,  and  22  to  25,  presents  the  results  of  six 
separate  computations  of  the  lift  ratio  as  a  function  of  the 

distance  penetrated  into  a  sharp-edge  gust  and  shows  considerable 
scatter.  The  derivation  given  by  Kiissner  in  reference  5  is  probably  the 
best.  No  direct  experimental  verification  of  this  factor  is  available, 
but  Kiissner  has  made  several  attempts  at  verification  by  dropping  air¬ 
foils  equipped  with  end  plates  through  the  boundary  of  a  wind-tunnel  jet. 
The  tests  made  by  Kiissner  have  been  described  in  reference  5  and  the 
comparisons  of  the  computed  and  experimental  flight-path  curvatures 
indicate  good  agreement,  at  least  for  penetrations  up  to  3  to  4  .chords. 
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The  values  of  Cp,  and  Cl  for  the  infinite-aspect -ratio  wing  are 
a  g 

presented  in  figure  26.  If  the  pitching  motion  of  the  airplane  is  not 
considered,  the  exact  shapes  of  the  unsteady-lift  functions  are  unimportant 
and  the  discrepancies  noted  in  figure  25  do  not  appear  to  he  significant. 

Finite  aspect  ratio. -  The  effect  of  finite  aspect  ratio  on  the 
unsteady-lift  functions  has  "been  investigated  ty  Jones  in  reference  14 
and  these  results  are  shown  in  figure  26.  This  figure  indicates  that, 
as  the  aspect  ratio  decreases,  the  lift  develops  more  rapidly,  parti¬ 
cularly  at  the  start,  until  for  the  limiting  case  of  zero  aspect  ratio 
the  variation  in  lift  might  he  expected  to  correspond  to  that  for  steady 
flow  conditions.  The  estimation  of  the  unsteady -lift  functions  for 
lower  aspect  ratios  than  those  shown  in  figure  26  has  teen  made  hy  extra¬ 
polating  the  data  ty  assuming  that  at  zero  aspect  ratio,  the  unsteady- 
lift  function  disappears.  Unpublished  tests  made  with  a  flying  wing 
having  an  aspect  ratio  of  about  I.27  showed  that  this  method  led  to 
reasonable  values.  The  calculated  value  of  the  acceleration  ratio 
was  O.92  as  compared  to  an  average  experimental  value  of  0*93*  These 
results  indicate,  therefore,  that  for  flying  wings  finite  theory  applies. 
Other  results,  shown  in  figure  27  and  given  by  Keuthe  in  reference  17 
and  by  Sears  and  Keuthe  in  reference  26,  indicate  that  the  unsteady  lift 
developed  on  a  low-aspect -ratio  wing  subjected  to  a  gust  more  closely 
approximates  the  unsteady-lift  function  for  the  two-dimensional  case 
than  that  for  aspect  ratio  3* 

An  indirect  verification  of  unsteady-lift  functions  has  been  obtained 
by  comparing  the  calculated  and  actual  response  of  airplane  models  to  a 
sharp  gust.  The  curves  in  figure  28  show  the  calculated  acceleration 
ratio  An/Ans  as  a  function  of  the  mass  parameter  Pg.  The  calculation 
for  aspect  ratio  6  is  based  on  the  unsteady -lift  function  given  in 
figure  26  and  that  for  infinite  aspect  ratio  is  based  on  the  unsteady- 
lift  functions  given  by  Rhode  in  reference  4.  The  test  points  shown  in 
figure  28  represent  experimental  values  of  the  acceleration  ratio  deter¬ 
mined  from  tests  in  both  the  old  and  the  new  gust  tunnels.  The  results 
are  for  all  types  of  airplanes  from  tailless  high-speed  configurations 
with  wings  swept  back  35°  to  transports.  As  can  be  seen  from  figure  28, 
the  experimental  data  scatter  about  the  curve  based  on  the  unsteady-lift 
functions  for  aspect  ratio  00  and  below  the  curve  for  aspect  ratio  6. 

The  results  indicate  that,  for  airplanes  that  include  a  fuselage,  the 
use  of  the  unsteady  functions  for  infinite  aspect  ratio  have  given 
adequate  and  in  mdst  cases  more  accurate  predictions  of  the  acceleration 
ratio  than  the  computations  based  on  the  finite-aspect -ratio  theory  of 
Jones. 

The  discrepancy  between  experiment  and  theory  was  considered  as  to 
the  effect  of  interference  on  the  experimental  results.  If  the  fuselage 
is  considered  as  an  elongated  wing  having  a  chord  three  times  that 
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cf  the  wing,  then  at  peak  acceleration  in  a  sharp  gust,  the  wing  has 
traveled  about  4  chords  while  the  fuselage  has  traveled  1.  On  this  basis, 
the  lift  on  the  fuselage  might  be  considered  zero  or  at  most  about  half 
of  its  steady-flow  value.  If  this  hypothesis  is  compared  with  the  mater¬ 
ial  given  in  reference  19  (some  of  the  values  were  in  error  and  have  been 
recomputed) ,  the  results  shown  in  table  VI  are  obtained.  The  results  in 
table  VI  indicate  that,  although  Jones'  unsteady  functions  are  much  too 
high  on  the  basis  of  gross  area  and  the  functions  from,  reference  4  are 
somewhat  high,  the  use  of  a  net  wing  area  plus  one-half  the  fuselage 
intercept  brings  all  calculations  into  closer  agreement.  In  fact,  for 
the  infinite-aspect -ratio  functions,  the  discrepancies  are  less  than  the 
precision  of  the  data.  If  the  net  wing  area  is  assumed,  then  both  finite- 
and  infinite-aspect -ratio  functions  differ  from  experiment  by  about  the 
same  amount,  the  results  based  on  reference  14  being  high  and  those 
based  on  reference  4  being  low*  Similar  corrections  may  apply  to  Keuthe's 
results  (fig.  27) .  Although  the  evidence  indicates  that  the  unsteady- 
lift  functions  for  infinite  aspect  ratio  should  be  used  for  conventional 
airplanes  in  conjunction  with  the  net  wing  area  plus  half  the  fuselage 
intercept,  the  use  of  net  area  is  recommended  and  is  discussed 
subsequently. 

Unpublished  tests  of  a  skeleton  airplane  equipped  Vith  a  wing  swept 
back  45°  showed  that  the  installation  of  the  fuselage  had  no  appreciable 
effect  on  the  maximum  acceleration  increment.  In  this  particular  case, 
however,  the  length  of  the  wing  from  the  leading  edge  of  the  root  section 
to  the  trailing  edge  of  the  tip  section  wa3  almost  equal  to  the  length 
of  the  fuselage.  Since  one  effect  of  sweep  would  be  to  modify  the  rate 
of  development  of  lift  on  the  wing,  the  difference  between  the  lift 
developed  on  the  wing  and  that  on  the  fuselage  for  a  given  gust  pene¬ 
tration  may  have  been  too  small  to  be  noted  during  the  test-  Although 
the  evidence  is  still  scant  and  conflicting  in  some  respects,  the  use  of 
net  wing  area  appears  to  be  better  than  the  use  of  gross  wing  area  for 
the  sharp-edge  gust  except  when  the  side  projections  of  the  wing  and  the. 
fuselage  are  about  the  same  length. 

The  determination  of  the  proper  wing  area  for  the  gradient-  gust  is 
complicated  by  the  introduction  of  the  pitching  motion  of  the  airplane, 
and  comparisons  must  be  made  between  the  detailed  calculations  and  experi¬ 
ment.  In  the  investigations  that  have  been  made  (reference  15)  the 
results  have  indicated  that  the  use  of  the  net  wing  area  yields  the  best 
over-all  agreement  between  calculations  and  experiment  for  gradient 
distances  between  0  and  1 6  chords. 

The  use  of  tapered  and  swept  wings  for  modern  aircraft  leads  to 
problems  in  the  computation  of  the  gust  load  factor.  Relatively  little 
information  is  available  on  either  problem,  but  the  data  of  reference  19 
indicated  that  moderate  amounts  of  taper  (up  to  2:1)  have  little  or  no 
effect  and  can  be  neglected  at  least  until  experimental  evidence  of  a 
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more  accurate  nature  is  available •  Large  values  of  taper  (for  example, 
of  about  3:l)  and.  large  sweep  angles  (30°  to  40°)  lead,  to  considerable 
concern  as  to  the  adequacy  of  the  unsteady-lift  functions,  particularly 
the  Kiissner  function  CT  •  The  concern  arises  from  the  fact  that  for 

wings  of  high  taper,  the  root  section  with  the  larger  chord  will  have  a 
rate  of  development  of  lift  for  a  given  distance  penetrated  into  a  gust 
considerably  less  than  that  of  the  narrower  tip  chord.  In  the  case  of 
sweep,  the  root  of  a  sweptback  wing  penetrates  the  gust  first  and  the 
lift  may  develop  an  appreciable  value  before  the  tip  sections  ever  enter 
the  gust. 

With  regard  to  the  effect  of  wing  taper  on  the  unsteady-lift 
function,  no  theoretical  developments  are  available  to  permit  accurate 
computation,  and  calculations  have  been  made  by  utilizing  strip  theory 
to  develop  the  unsteady-lift  functions  for  the  finite  wing  from  the  two- 
dimensional  functions  of  Kiissner  and  Wagner.  Experimental  data  available 
to  check  the  effect  of  taper  on  the  unsteady-lift  functions  are  indirect 
and  are  the  result  of  tests  of  two  specific  airplanes.  In  both  cases, 
calculations  based  on  the  results  presented  in  reference  4  were  in  agree¬ 
ment  with  the  experimental  data  for  the  sharp -edge  gust  in  which  pitch 

can  be  neglected.  In  the  case  of  the  large  flying  boat  with  taper 

ratio,  the  discrepancy  amounted  to  about  2  percent.  The  maximum  dis¬ 
crepancy  was  well  within  the  experimental  error.  On  the  basis  of  these 

limited  results,  the  effect  of  taper,  at  least  up  to  3^:1,  appears  to  be 
negligible  insofar  as  the  calculation  of  total  loads  is  concerned. 

Wo  theoretical  studies  are  available  for  the  sweptback  wing,  but 
recent  test  results  for  a  straight  and  a  45°  sweptback  wing  are  available. 
Tests  were  made  on  a  straight  wing  with  a  2:1  taper  and  on  the  equivalent 
of  the  straight  wing  where  each  half-wing  was  rotated  about  the  midchord 
point  at  the  root  so  that  the  span  changed  with  the  angle  of  sweep.  The 
sweep  was  such  that  the  midchord  line  was  at  an  angle  of  45°  to  its 
original  position.  Flights  were  made  through  a  sharp-edge  gust,  and  an 
average  time  history  of  acceleration  increments  as  a  fraction  of  the 
maximum  value  is  shown  in  figure  29.  For  comparison  with  the  experi¬ 
mental  results,  two  curves  are  shown,  one  based  on  the  theory  presented 
in  reference  4  which  disregards  the  effect  of  sweep  on  the  unsteady -lift 
functions  and  the  other  represents  the  calculated  time  history  of  accele¬ 
ration  based  on  the  assumption  that  strip  theory  could  be  applied  to  derive 
the  Kiissner  function  for  the  swept  wing. 

Figure  29  shows  that  strip  theory  is  in  good  agreement  with  experi¬ 
ment  throughout  the  entire  history,  but  the  neglect  of  the  effect  of 
sweep  on  the  lag  in  lift  is  in  error  at  the  start  of  the  motion,  although 
good  agreement  is  obtained  for  distances  greater  than  4  chords.  Both 
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calculations  should  yield  essentially  the  same  values  for  the  larger  . 
distances  since  the  conditions  tend  to  approach  the  steady  state.  The 
discrepancy  between  experiment  and  strip  theory  appears  to  be  due  in  part 
to  positive  pitch  of  the  model  during  the  traverse  of  the  gust.  Therefore, 
the  effect  of  sweep  on  the  unsteady-lift  functions  cannot  be  neglected 
if  the  shape  of  the  curve  is  considered  important. 

The  unsteady-lift  theory  as  developed  by  Wagner,  Jones,  and  others 
has  been  evolved  for  monoplanes  only  but  the  unsteady -lift  functions  for 
biplanes  may  be  needed  in  connection  with  gust -load  calculations.  Wo 
theory  is  available  on  this  problem,  but  consideration  of  the  vortex 
sheets  associated  with  the  biplane  indicates  that,  for  -the  first  few 
chords  after  a  sudden  change  in  angle  of  attack,  the  tip  vortices  of  the 
wings  are  short  and  have  a  negligible  effect  on  the  mutually  induced 
angle  of  attack  of  one  wing  upon  the  other.  Since  the  equivalent  mono¬ 
plane  theory  assumes  completely  developed  vortex  systems  with  the  shed 
vortex  at  infinity,  it  might  be  expected  that  for  a  biplane  in  a  sharp 
gust  the  two  wings  would  act  independently,  and,  therefore,  the  unsteady- 
lift  functions  and  other  associated  aerodynamic  parameters  for  the  bi¬ 
plane  should  be  based  on  the  characteristics  for  each  individual  wing. 
Figure  30,  which  has  been  reproduced  from  reference  8,  shows  the  results 
of  tests  of  a  biplane.  The  results  in  the  figure  are  the  accelerations 
obtained  as  a  result  of  experiment  in  the  Langley  gust  tunnel,  the  curve 
predicted  by  assuming  that  the  two  wings  act  independently,  and  the  curve 
predicted  by  assuming  that  the  equivalent-monoplane  theory  holds. 
Infinite-aspect -ratio  unsteady-lift  functions  were  used.  For  the  shorter 
gust-gradient  distances  where  the  tip  vortices  might  be  expected  to  have 
little  effect,  theory  and  experiment  are  in  excellent  agreement  if  the 
two  wings  are  assumed  to  act  independently.  The  computed  curve  for 
acceleration  increment,  based  on  reference  k,  is  considerably  below  the 
experimental  dataj  thus  the  tip  vortices  have  little  effect,  at  least 
up  to  8  or  10  wing  chords.  As  the  gradient  distance  is  increased  to 
20  chords,  the  experimental  data  fall  below  both  calculated  curves  because 
of  the  pitching  action  of  the  airplane. 

The  effect  of  compressibility  on  the  unsteady -lift  functions  is  of 
interest  because  it  may  change  the  shape  of  the  curves,  but  little  or 
no  information  on  this  problem  is  available.  Jones  has  indicated  that 
for  subsonic  speeds  the  effect  of  compressibility  would  be  to  reduce  the 

aspect  ratio  of  the  wing  according  to  the  relation  \j  1  -  M^.  The  E 
correction  of  reference  27  is  changed  at  the  same  time  to  correspond. 

For  the  transonic  range  where  mixed  flows  occur  the  variation  in  the 
unsteady-lift  functions  is  unknown.  At  supersonic  speeds  the  effect 
of  unsteady  lift  would  be  expected  to  disappear  since  the  effect  requires 
the  transmission  of  pressures  forward.  It  would  be  expected  that  the 
unsteady-lift  functions  would  change  only  slightly  up  to  high  subsonic 
speeds  and  that  at  supersonic  speeds,  the  functions  would  disappear. 
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Slop©  of  Lift  Curve 

Investigations  to  determine  the  proper  slope  of  the  lift  curve  to 
he  used  in  gust-load  calculations  have  been  made  mainly  in  the  Langley 
gust  tunnel.  The  information  available  for  presentation  is  indirect 
and  in  many  instances  inconclusive  because  of  interrelations  of  many 
factors.  The  factors  involved  in  the  computation  of  the  proper  slope  of 
the  lift  curve  include  the  determination  of  the  airfoil  section  character¬ 
istics,  aspect-ratio  corrections,  interference  effects,  the  effects  of 
configuration  with  reference  to  sweep  and  multiplanes  and  the  effects  of 
compressibility  and  power. 

Section  characteristics. -  As  indicated  by  the  results  presented  in 
figure  28,  the  conventional  estimates  of  the  lift-curve  Blope  together 
with  the  unsteady -lift  functions  for  infinite  aspect  ratio  gave  the  best 
results  in  predicting  the  acceleration  increment  due  to  a  sharp-edge 
gust  for  ordinary  airfoils.  The  use  of  the  laminar-flow  section  intro¬ 
duced  airfoils  having  higher  slopes  of  the  lift  curve.  The  increase 
was  of  considerable  concern  since  it  amounted  to  an  increase  of  10 
percent  in  the  load-factor  increment.  Brief  consideration  indicated 
that  the  slope  of  the  lift  curve  depends  on  the  steady-flow  relation 
between  the  boundary-layer  thickness  and  angle  of  attack.  It  was 
believed  that  under  unsteady  flow  conditions  the  relation  for  steady 
flow  would  not  apply  and,  therefore,  the  boundary  layer  would  not  have 
time  to  adjust  itself  during  a  sudden  change  in  angle  of  attack.  Since 
no  experimental  or  theoretical  information  was  available  in  connection 
with  the  probable  slope  of  the  lift  curve,  experimental  evidence  was 
needed. 

For  this  purpose  a  skeleton  model  with  a  laminer-flow  wing  was 
tested  in  the  Langley  gust  tunnel.  The  characteristics  of  the  test 
model  are  given  in  table  VII  and  figure  31*  For  the  first  tests  the 
wing  had  a  smooth  surface.  For  the  second  tests,  carborundum  grains 
were  glued  over  the  leading  edge  of  the  wing  back  to  7*8  percent  of  the 
chord  on  both  the  upper  and  lower  surfaces.  The  slopes  of  the  lift 
curve  for  the  smooth  and  rough  wings  have  been  included  in  table  VII  and 
are  based  on  experimental  data  presented  in  reference  28  corrected  to 
finite  aspect  ratio  according  to  reference  27.  In  addition  to  the  lift- 
curve  slopes  obtained  from  low-turbulence  wind-tunnel  tests,  the  theore¬ 
tical  lift-curve  slope,  based  on  the  theory  presented  in  reference  29,  is 
also  included  in  the  table. 

The  average  values  of  the  acceleration  increments  for  the  smooth 
and  rough  wings  and  the  probable  errors  are  given  in  table  VIII.  The 
acceleration  increments  computed  on  the  basis  of  the  characteristics  in 
table  VII  have  also  been  included.  Table  VIII  shows  that  the  maximum 
acceleration  increments  for  the  two  test  conditions  are  essentially  the 
same.  Comparison  of  the  0.02g  experimental  difference  with  the  calculated 
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difference  of  0.26g  shows  that  roughening  an  airfoil  to  produce  a 
turbulent  boundary  layer  has  no  effect  on  the  slope  of  the  lift  curve 
that  is  applicable  in  the  unsteady  flow  conditions  of  a  sharp-edge  gust. 
If  these  results  are  extended  by  the  assumption  that  the  flow  conditions 
for  a  conventional  airfoil  in  the  steady  state  are  simulated  by  those 
for  the  roughened  low-drag  wing,  airfoil  section  characteristics  are 
concluded  to  have  no  significant  effect  on  the  acceleration  increment 
obtained  in  unsteady  flow  conditions  of  a  sharp-edge  gust. 

Since  the  most  probable  gust  is  one  with  a  gradient  distance  of 
10  chords,  the  averages  of  the  maximum  acceleration  increments  for  the 
two  test  conditions  were  calculated  for  gusts  with  gradient  distance 
up  to  12  chords  by  applying  the  principle  of  superposition  to  the  test 
results  for  the  sharp-edge  gust.  The  results  of  these  calculations 

showed  only  2^- -percent  difference  between  the  average  values.  On  the 

basis  of  this  simple  analysis,  airfoil  section  characteristics  have  no 
significant  effect  on  the  acceleration  increments  for  gusts  with  grad¬ 
ient  distances  from  0  to  12  chords.  From  the  present  data,  however,  it 
is  not  possible  to  specify  the  gradient  distance  at  which  the  difference 
would  become  noticeable  an  the  result  of  attaining  quasi -steady 
conditions. . 

For  the  laminar-flow  type  of  airfoil,  tests  in  low-turbulence  wind 
tunnels  together  with  the  aspect -ratio  corrections  of  reference  27 
should  yield  an  adequate  prediction  of  the  acceleration  increment  due  to 
a  sharp-edge  gust.  The  use  of  the  corrections  and  the  determination  of 
section  characteristics  for  an  arbitrary  airfoil  are,  however,  still 
open  to  question  since  tests  in  low-turbulence  tunnels  on  the  arbitrary 
airfoil  sections  tend  to  show  low  slopes  of  the  lift  curve.  In  view  of 
the  many  interrelated  factors,  some  of  which  are  discussed  in  subsequent 
paragraphs,  it  appears  that  the  selection  of  an  arbitrary  section  lift- 
curve  slope  of  about  6  per  radian  will  be  adequate  for  gust -load 
calculations. 


Aspect-ratio  corrections.  -  As  previously  noted,  the  data  in  fig¬ 
ure  28  were  obtained  for  conventional  airplanes,  but  the  data  in  table  VIII 
apply  to  a  skeleton  airplane  with  little  or  no  fuselage.  The  combined 
effects  of  fuselage  interference,  aspect-ratio  corrections,  unsteady- 
lift  functions,  and  section  characteristics  are  all  involved  and  at  the 


present  time  cannot  be  entirely  segregated.  For  the  time  being,  simple 


aspect -ratio  correction  factors  appear  to  be  adequate. 


The  factor 


A  +  2 


has  been  used  with  satisfactory  results  in  connection  with  gust-load 
studies. 
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Swept  wings . -  Table  IX  (from,  reference  30)  shows  the  results  of 
gust-tunnel  tests  on  a  wing  swept  back  45°.  The  tests  were  made  in  a 
sharp-edge  gust  and  the  acceleration  increments  for  a  gradient  distance 
of  9  chords  were  obtained  by  superposition-  All  results  are  corrected 
to  zero  pitch  of  the  model.  Calculated  values  are  also  included.  One 
set  of  values  is  based  on  the  results  of  wind-tunnel  tests  of  the  air¬ 
plane  model,  tail  off,  and  the  other  set  differs  from  the  first  in  that 
the  slope  of  the  lift  curve  for  the  swept  wing  wan  taken  as  that  for  the 
straight  wing  multiplied  by  the  cosine  of  the  sweep  angle.  In  both  cal¬ 
culations,  the  unsteady -lift  function  for  the  wing  was  obtained  by 

means  of  strip  theory.  8 

Inspection  of  the  results  given  in  table  IX  shows  that  the  use  of 
the  cosine  law  for  predicting  the  slope  of  the  lift  curve  for  the  swept 
wing  gives  excellent  results,  the  differences  between  experiment  and 
calculation  being  within  the  experimental  error.  The  use  of  lift -curve 
slopes  from  steady-flow  tests  to  calculate  the  acceleration  increments 
for  the  swept  wing  yields  values  some  20  percent  below  experimental 
results.  Similar  tests  on  a  4-5°  sweptforward  wing  verify  this  result. 

The  low  value  of  acceleration  increments  calculated  from  the  results 
of  wind-tunnel  tests  has  been  of  concern.  At  present,  no  definite 
evidence  is  available  to  explain  this  discrepancy,  but  it  is  believed 
that  the  difference  can  be  ascribed  to  the  behavior  of  the  boundary 
layer  in  the  unsteady  flow  condition  that  exists  during  traverse  of  a 
gust.  The  discrepancy  emphasizes  the  importance  in  the  prediction  of 
gust  load  factors  of  using  the  proper  lift-curve  slope.  On  the  basis 
of  limited  data  better  results  seem  to  be  obtained  if  the  cosine  law  of 
variation  is  used  rather  than  wind-tunnel  results. 

The  results  obtained  to  date  apply  to  a  definite  method  of  sweeping 
the  wing  (rotating  the  wing  panel)  to  obtain  the  equivalent  straight 
wing.  Other  systems  lead  to  other  combinations  of  aspect  ratio  and  sweep 
and  may  yield  different  results.  In  such  calculations,  the  same  procedure 
as  that  described  in  the  present  paper  should  be  utilized. 

Scale  effects. -  Tests  in  the  old  and  new  gust  tunnels  of  models  of 
the  Boeing  B-247  airplane  (table  X)  were  made  to  obtain  a  measure  of  any 
scale  effects.  A  3"foot-span  model  was  tested  in  the  old  gust  tunnel  as  a 
reference  for  design  requirements.  When  the  new  gust  tunnel  was  built 
this  model  and  one  dynamically  similar  but  twice  as  large  were  tested. 

The  results  of  all  three  tests  are  shown  in  figure  32  as  a  plot  of  the 
acceleration  ratio  as  a  function  of  gradient  distance  in  chords. 

Inspection  of  figure  32  shows  excellent  agreement  between  the 
results  for  the  two  models^  this  fact  indicates  that  for  the  conditions 
tested  no  scale  effect  was  present.  The  largest  discrepancy  is  for  a 
sharp-edge  gust  and  is  no  greater  than  differences  in  results  for  tests 
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of  the  same  model  in  two  tunnels.  Same  of  the  discrepancy  may  he  due  to 
the  fact  that  the  gust  shape  was  the  same  in  absolute  dimension  hut 
differed  slightly  on  the  basis  of  chords  of  the  two  models. 

Effect  of  -power. -  When  the  slipstream  covers  the  entire  span  of  the 
airplane,  the  steady-flow  slope  of  the  lift  curve  has  been  shown  to 
increase  about  100  percent  for  the  slipstream  corresponding  to  the  climb 
condition.  Since  the  working  velocity  has  been  increased  by  the  intro¬ 
duction  of  power,  a  real  increase  in  the  lift-curve  slope  has  probably 
been  obtained.  Under  these  circumstances,  the  angle -of -attack  change  due 
to  the  gust  would  be  utilized  in  connection  with  the  increased  slope  of 
the  lift  curve.  In  actual  practice,  the  suggestion  that  the  power-on 
slope  of  the  lift  curve  should  be  used  where  the  slipstream  covers  most 
of  the  airplane  span  is  not  as  serious  as  it  appears  at  first  sight. 

The  highest  lift-curve  slope  is  obtained  at  relatively  low  forward  speed, 
a  condition  at  which  the  load  factor  due  to  the  gust  is  small.  For 
conventional  airplanes,  the  power-off  lift-curve  slope  appears  correct. 

Multiplanes. -  For  steady  flow  conditions  the  biplane  is  usually 
represented  by  an  equivalent  monoplane.  This  assumption  no  longer 
applies  for  the  transient  lift  conditions  that  exist  in  a  gust.  The 
question  of  the  proper  slope  of  the  lift  curve  to  use  for  the  transient 
conditions  is  in  part  answered  by  the  results  presented  In  figure  30. 

On  the  basis  of  the  agreement  shown  in  figure  30,  the  wings  should  be 
considered  as  acting  independently  and  the  slope  of  the  wing  lift  curve 
as  being  computed  on  the  basis  of  the  average  geometric  aspect  ratio  of 
the  two  wings. 

Compressibility ■ -  The  effect  of  compressibility  on  the  slope  of  the 
lift  curve  under  unsteady  flow  conditions  is  a  subject  of  much  Interest 
and  one  for  which  no  experimental  evidence  Is  available-  It  is  thought 
that  the  lift -curve  slope  should  follow  the  Glauert  factor  in  the 
subsonic  range,  but  in  the  region  near  the  critical  Mach  number  and 
through  the  transonic  region  the  mixed  flows  combined  with  transient 
conditions  might  delay  or  cancel  the  effect  of  comperes sibility  on  the 
lift-curve  slope.  For  the  time  being  the  usual  compressibility  factors 
could  be  applied  up  to  the  critical  Mach  number  and  above  that,  the 
results  of  wind-tunnel  tests  should  be  utilized. 


Downwash 

Figure  33  (fig-  9  of  reference  31)  shows  the  change  in  vertical 
velocity  at  the  tail  following  a  sudden  increase  in  wing  circulation  for 
several  values  of  the  distance  from  the  trailing  edge  of  the  wing  to  the 
leading  edge  of  the  tail.  The  change  in  vertical  velocity  is  plotted  as 
a  function  of  the  distance  of  the  tail  surface  from  the  vortex  shed  by 
the  wing. 
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Figure  33  indicates  that  the  tail  surface  first  encounters  a  grad¬ 
ual  increase  in  upwash  until  the  leading  edge  crosses  the  center  of  the 
shed  vortex,  vhen,  following  a  violent  change  in  direction,  the  downwash 
approaches  its  steady-state  value.  The  curves  for  different  tail  lengths 
indicate  quite  clearly  that  the  effect  of  increasing  the  tail  length  is 
to  increase  the  time  delay  of  the  velocity  change  in  proportion  to  the 
tail  length.  The  violence  of  the  direction  change  as  the  tail  penetrates 
the  vortex  depends  on  the  vertical  location  of  the  tail  relative  to  the 
vortex,  and,  until  more  accurate  predictions  can  he  made  of  the  location 
of  the  vortex,  such  changes  might  he  disregarded.  Figure  33  indicates 
that  the  assumption  proposed  in  references  31  and.  32  of  a  space  lag 
between  the  change  in  lift  on  the  wing  and  the  change  in  angle  of  attach 
at  the  tail  is  reasonable  and  should  he  taken  into  account  in  any 
detailed  calculations. 

The  indications  to  date  are  that  the  results  obtained  by  utilizing 
a  space  lag  of  downwash  have  been  in  fair  agreement  with  experiment,  but 
the  corresponding  calculations  neglecting  the  space  lag  of  downwash  have 
not  been  made  to  determine  the  seriousness  of  the  error  of  neglecting 
this  quantity. 


Maximum  Lift  Coefficient 

Reference  19  shows  that  for  an  airplane  model  traversing  a  sharp- 
edge  gust  at  its  steady-flcw  maximum  lift  coefficient,  the  maximum  lift 
coefficient  during  the  action  of  the  gust  was  not  limited  to  the  steady- 
flow  value*  Farren  in  reference  33  made  tests  of  two-dimensional  air¬ 
foils  at  constant  angular  velocities  through  maximum  lift  and  return. 

Angle-of-attack  variations  range  as  high  as  12^  per  chord  of  travel  and 

the  results  indicated  that  the  maximum  lift  coefficient  could  range  from 
30  to  50  percent  above  the  steady-flow  value. 

No  concise  estimates  of  maximum  lift  coefficient  during  sudden 
changes  in  angle  of  attack  can  be  made,  but  the  sketchy  information  that 
is  available  indicates  that  a  conservative  estimate  for  incompressible- 
flow  conditions  is  an  increase  of  about  25  percent  over  the  steady-flow 
value.  In  the  transonic  region,  the  effects  of  compressibility  may 
limit  the  maximum  lift  coefficient,  and  in  this  range  the  investigations 
made  in  the  various  wind  tunnels  would  be  pertinent. 


RIGID-BOUT  REACTIONS 


The  available  data  op  gust  structure  and  aerodynamics  have  been 
analyzed  and  the  next  problem  is  the  determination  of  the  behavior  of 
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the  airplane  as  a  whole  since  the  motions  of  the  airplane  have  an 
Important  hearing  on  gust -load  calculation.  The  significant  parameters 
of  airplane  behavior  as  determined  from  analytical  and  experimental 
investigations  are  subsequently  summarized  and  possible  simplifications 
investigated. 


Analytical  and  Experimental  Studies 

The  material  presented  is  the  result  of  analytical  and  experimental 
study  of  arbitrary  and  special  configurations.  The  analytical  studies 
have  been  of  a  general  nature  and  about  two-thirds  of  the  experimental 
studies  can  be  so  classed.  The  rest  of  the  experimental  research  has 
been  done  in  connection  with  specific  problems  or  airplane  designs.  The 
scope  and  results  of  the  analytical  studies  are  presented  first,  followed 
by  the  experimental  studies.  In  addition  to  the  behavior  of  and  the 
loads  on  airplanes  caused  by  vertical  gusts  uniform  across  the  span, 
limited  research  and  studies  on  unsymmetrical  and  lateral  gusts  are 
described  and  reported  in  appropriate  sections. 

Analytical  studies. -  The  analytical  studies  of  the  loads  on  and 
behavior  of  an  airplane  traversing  a  gust  have  been  performed  for  a 
uniform  upward-acting  vertical  gust.  Current  theory  indicates  that  the 
incremental  loadings  for  a  downward-acting  gust  are  equal  in  magnitude 
but  opposite  in  sign.  Further  assumptions  are  that: 

(1)  The  airplane  maintains  a  constant  forward  speed  during  the 
traverse  of  the  gust 

(2)  The  airplane  is  in  equilibrium  prior  to  entry  into  the  gust 

(3)  The  aerodynamic  center  of  an  aerodynamic  surface  is  at  the 
quarter-chord  point 

(4)  The  moment  coefficient  at  zero  lift  is  a  constant  for  transient 
conditions  and  has  the  same  value  as  for  steady  flow  conditions 

(5)  The  control  surfaces  of  the  airplane  are  locked 

Extended  analyses  in  which  the  airplane  is  considered  free  to  pitch 
as  well  as  to  rise  under  the  action  of  a  gust  have  been  made  by  the  NACA 
for  conventional,  canard,  and  tailless  airplanes.  In  the  case  of  the 
conventional  airplane,  the  analysis  was  rather  wide  in  scope;  whereas 
for  the  other  two  types  of  aircraft,  the  analyses  were -limited  to  specific 
cases.  The  procedure  used  was  an  iteration  process. 

The  calculations  were  aimed  at  covering  all  configurations  and  values 
of  stability  which  would  be  reasonable  for  conventional  airplanes  and 
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such  values  of  the  mass  parameter  and  moment  of  inertia  vhich  would  fall 
within  the  capacity  of  the  gust  tunnel.  The  pertinent  airplane  charac¬ 
teristics  assumed  in  the  analysis  are  given  in  table  XI  for  all  the  com¬ 
binations  considered.  As  can  be  seen  from  the  table,  the  static  margins 
range  from  0.009  to  -0-503  and  the  mass  parameter  ranged  from  10. 25  to 
30-75*  The  calculations  were  made  for  three  gust-gradient  distances  - 
0,  8,  and  l6  chords.  The  flight  conditions  assumed  were  a  gust  velocity 
of  6  feet  per  second  and  a  forward  speed  of  about  40  miles  per  hour.  The 
calculations  were  not  made  for  all  combinations  of  variables  noted  in 
table  XI,  but  same  intermediate  values  were  interpolated. 

The  results  of  the  calculations  are  shown  in  tables  Xll(a),  Xll(b), 
and  XIl(c)  as  the  acceleration  increments  at  maximum  acceleration  for 
the  airplane  with  the  center  of  gravity  at  15,  25,  and  35  percent  of  the 
mean  aerodynamic  chord,  respectively.  Each  table  shows  the  results  for 
the  three  gust-gradient  distances  and  for  various  combinations  of  tail 
area  and  tall  length.  The  incremental  values  of  acceleration  which  are 
contributed  by  each  motion  or  source  are  tabulated  and  then  totaled  to 
obtain  the  wing  acceleration,  the  acceleration  Increment  resulting  from 
the  lift  on  the  tail,  and  finally  the  total  acceleration  increment 
impressed  on  the  airplane. 

In  addition  to  the  results  of  the  analysis  being  presented  in 
tabular  form,  figures  3^  and  35  indicate  the  effect  of  pitch  on  the  total 
wing  load.  In  figure  3^>  the  ratio  of  the  total  wing  load  increment, 
including  the  effect  of  pitch,  divided  by  the  total  wing  load  increment, 
assuming  the  pitch  equal  to  zero,  is  shown.  A  value  of  1  indicates  that 
the  effect  of  pitch  was  negligible  and  values  greater  than  1  indicate 
that  the  effect  of  pitch  Increased  the  total  wing  load.  The  results  are 
shown  as  a  function  of  the  static  margin  for  three  gust-gradient  dis¬ 
tances.  The  individual  curves  shown  In  the  figure  represent  given  con¬ 
figurations  but  different  center-of -gravity  positions.  Figure  35  is  a 
similar  plot  of  data  from  tables  Xll(a),  Xll(b),  and  XIl(c)  of  the 
total  acceleration  increment  as  a  function  of  the  static  margin  for 
comparison  with  the  experimental  data  from  the  gust  tunnel. 

Figure  36  shows  the  wing  lift  increment  resulting  from  pitch  of  the 
airplane  In  terms  of  the  total  wing  lift  as  a  function  of  mass  parameter. 
The  three  solid  curves  represent  the  three  center -of -gravity  positions. 
The  curves  represent  the  "average"  airplane  with  medium  tail  length  and 
medium  tail  area  traversing  a  flat-top  gust  with  a  gradient  distance  of 
8  chords.  The  dash  lines  included  in  the  figure  represent  the  pitch- 

A0 

increment  ratio  — r*.  This  ratio  has  been  used  to  correct  for  pitch 

U/V 

effects  on  the  basis  that  the  change  In  acceleration  due  to  pitch  is 
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proportional  to  the  ratio  of  the  pitch  increment  to  the  gust  angle.  The 
distance  between  the  dash  and  solid  lines  is  a  measure  of  the  error  of 
this  assumption. 

The  data  in  tables  Xll(a),  Xll(b),  and  XIl(c)  were  also  utilized  to 
obtain  the  total  tail  load  increment  divided  by  the  tail  load  if  the 
downwash,  the  vertical  motion,  and  the  pitch  are  assumed  to  be  zero. 

The  results  are  shown  in  figure  37  as  a  function  of  the  static  mar¬ 
gin  dCm  /  dCr  for  the  three  gust -gradient  distances.  The  total  tail 
eg/ 

load  shown  is  that  at  the  time  of  maximum  total  airplane  load.-  The  lines 
connecting  the  points  represent  the  same  configuration  and  the  movement 
along  the  line  to  the  right  indicates  increasing  static  stability. 

Detailed  calculations  have  also  been  made  for  a  flying  wing  with 
about  25°  of  sweepback  and  for  a  canard  airplane-  The  characteristics 
of  the  airplanes  are  shown  in  table  XI.  The  calculations  for  the  flying 
wing  were  made  for  two  center-of -gravity  positions  and  for  three  gradient 
distances  -  0,  8,  and  17*5  chords-  Table  XIII  shows  the  results  of  the 
calculations  and  includes  the  contributions  from  the  different  sources 
for  the  two  center-of -gravity  positions.  The  total  wing  load  increments 
are  shown  for  each  center-of -gravity  position.  Calculations  were  made 
by  using  finite-aspect -ratio  unsteady-lift  functions  on  the  basis  that 
no  fuselage  was  present-  The  total  acceleration  increment  has  been 
plotted  in  figure  38(a)  for  each  center-of -gravity  position.  The  canard 
airplane  had  the  same  general  characteristics  as  the  conventional 
Boeing  B-2k7  transport  airplane-  The  calculations  for  the  canard  air¬ 
plane  for  net  wing  area  and  for  the  two  unsteady-lift  functions  have  been 
obtained  from  reference  15  and  are  presented  in  table  XIV .  The  calculated 
total  acceleration  increments  for  the  canard  are  shown  in  figure  38(b) 
as  a  function  of  gradient  distance. 

Figure  39  indicates  the  effect  on  the  calculated  loads  of  substitut¬ 
ing  the  wing  loading  for  the  mass  parameter.  This  figure  is  a  plot  of 
the  acceleration  based  on  wing  loading  divided  by  that  based  on  the  mass 
parameter  as  a  function  of  wing  loading.  The  calculations  were  made  for 
three  classes  of  airplanes  -  transports,  personal  airplanes,  and  flying 
boats.  The  results  are  shown  as  a  function  of  the  wing  loading.  Points 
lying  above  the  line  indicate  overestimation  of  the  acceleration  incre¬ 
ment,  and  points  below  the  line  represent  underestimation  of  the  accele¬ 
ration  increment. 

As  mentioned  in  the  section  about  gust  structure,  the  suggestion 
was  made  of  either  a  triangular  or  sinusoidal  gust  shape  in  the  direction 
of  flight.  As  a  check  on  the  interchangeability  of  the  two  gust  shapes, 
figure  kO  shows  the  ratio  of  the  acceleration  due  to  a  sinusoidal  gust 
to  that  for  a  triangular  gust  as  a  function  of  the  mass  parameter  when 
the  unsteady-lift  functions  for  infinite  aspect  ratio  and  aspect  ratio  6 
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are  used.  The  calculations  were  performed  by  a  step-by-step  procedure 
that  took  into  account  the  possible  lack  of  coincidence  in  the  gust  peaks 
Other  calculations  indicated  that  for  gusts  with  a  distance  of  10  chords, 
the  maximum  lag  in  the  acceleration  peak  would  be  about  2  chords  when 
the  wing  had  an  infinite  mass  parameter. 

Experimental  studies. -  Most  of  the  experimental  data  available  are 
from  general  tests  in  the  Langley  gust  tunnel  with  upward-acting  gusts 
perpendicular  to  the  flight  path.  Most  of  the  research  can  be  classed 
as  general  and  includes  the  material  presented  in  reference  19-  The 
results  of  tests  of  specific  models  that  are  applicable  to  general 
studies  are  included. 

A  general  study  was  made  of  the  effect  of  the  stability  of  airplanes 
on  the  gust  load  factor;  the  stability  characteristics  were  varied  over 
a  wide  range  within  which  conventional  airplanes  would  be  expected  to 
fall.  The  tests  were  performed  on  an  arbitrary  "stability"  airplane 
model  for  correlation  with  the  extended  analysis  previously  described. 

The  characteristics  of  this  model  are  given  in  table  XI  and  the  total 
acceleration  increments  obtained  from  the  tests  are  shown  in  figure  35* 
The  tests  were  made  for  three  gust -gradient  distances  and  with  the  center 
of-gravity  positions  at  17- 5>  25,  and  30  percent  of  the  mean  geometric 
chord.  The  combinations  of  tail  area  and  tail  length  used  on  the  test 
model  are  listed  in  figure  35. 

The  unconventional  airplane  models  listed  in  table  XI  include  a 
flying  wing  and  a  canard  airplane-  Tests  of  the  flying  wing  were  made 
with  the  center  of  gravity  at  20  percent  of  the  mean  geometric  chord. 

The  acceleration  increment,  forward  speed,  gust  velocity,  and  pitch  of 
the  model  were  measured  during  each  test  flight.  The  tests  consisted 
of  a  series  of  five  or  more  flights  of  the  airplane  at  a  forward  speed 
of  AO  miles  per  hour  through  vertical  gusts  with  gradient  distances  of  0} 
8,  and  16  chords.  The  test  results  have  been  included  in  table  XIII  and 
are  shown  in  figure  38(a).  The  test  conditions  for  the  canard  airplane 
have  previously  been  reported  in  reference  15-  The  pertinent  charac¬ 
teristics  of  the  model  have  been  included  in  table  XI  and  the  results  of 
the  experiments  are  given  in  table  XIV  and  in  figure  38(b) . 

A  few  flight  tests  have  been  made  to  obtain  data  on  the  behavior  of 
an  airplane  subjected  to  unsymmetrical  gusts  and  on  the  gust  loads  on 
the  vertical  and  horizontal  tail  surfaces.  The  results  for  the  unsym¬ 
metrical  gusts  have  been  discussed  previously  in  the  section  entitled 
"The  Structure  of  Atmospheric  Gusts  '  and  indicate  reasonable  agreement 
between  calculation  and  experiment.  The  results  of  investigations  of 
gust  tail  loads  on  the  0-2H  and  XB-15  airplanes  are  given  in  table  XV  as 
ratios  of  the  effective  gust  velocity  on  a  tail  surface  to  that  deter¬ 
mined  for  the  wing.  The  XB-15  airplane  Is  shown  in  figure  6  and  des¬ 
cribed  In  the  section  on  gust  structure  and  the  general  characteristics 
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of  the  0-2H  airplane  can  he  found  in  reference  3^*  Because  of  the  use 
of  improved  instruments,  the  results  from  work  with  the  XB-15  airplane 
are  considered  tetter  than  those  obtained  on  the  0-2H  airplane. 


Discussion 

Wing  area. -  The  results  presented  in  figures  35  and  38(B)  show  good 
agreement  between  experiment  and  the  detailed  calculation  based  on  net 
wing  area.  If  gross  area  were  used,  the  calculated  values  would  be 
higher  than  the  experimental  values  for  a  sharp  gust.  For  a  gradient 
gust  the  reverse  would  be  true  since  the  amount  of  alleviation  due  to 
pitch  would  be  increased.  A  similar  result  applies  to  the  canard  air¬ 
plane  of  figure  38(b)  as  indicated  by  table  II  of  reference  15,  which 
shows  a  decrease  in  acceleration  of  30  percent  for  a  gradient  distance 
of  16  chords  (A  =  °°  ) . 

Mass  parameter. -  Analytical  studies  that  are  based  on  extensive 
equations  indicate  that  the  mass  parameter  is  the  significant  variable 
in  that  it  affects  the  amount  of  alleviation  due  to  vertical  motion. 

The  mass-parameter  term  has  a  negligible  effect  for  a  sharp-edge  gust 
but  is  important  for  long  gradient  distances.  This  result  is  due  to  the 
fact  that  for  a  sharp  gust,  the  airplane  does  not  acquire  much  vertical 
velocity,  but  in  a  long  gust  the  vertical  velocity  approaches  the  gust 
velocity. 

In  addition  to  its  effect  on  the  vertical  motion  of  the  airplane, 
an  increase  in  the  mass  parameter  for  a  given  airplane,  by  increasing 
the  gross  weight,  unexpectedly  increases  the  alleviation  of  the  total 
wing  load  due  to  pitch.  (See  fig.  38.)  As  the  mass  parameter  increases, 
the  resistance  of  the  airplane  to  rotation  would  be  expected  to  increase 
and  therefore  the  amount  of  pitch  and  the  alleviation  due  to  the  pitching 
motion  would  be  decreased.  The  result  Just  noted  arises  from  the  fact 
that  the  pitching  moments  on  the  airplane  arise  not  only  from  the  action 
of  the  gust  but  also  from  the  vertical  motion  of  the  airplane.  By 
reducing  the  vertical  motion  of  the  airplane,  the  adverse  pitching  moment 
is  decreased,  and  thus  the  favorable  effect  predominates  and  provides 
more  alleviation. 

In  some  cases,  such  as  in  computing  the  factor  K,  wing  loading 
instead  of  mass  parameter  has  been  taken  as  the  significant  variable. 

Such  a  substitution  presumes  a  fixed  relation  between  wing  loading  and 
mass,  parameter.  Figure  39  indicates  that  this  substitution  may  result 
in  either  underestimation  or  over estimation  of  imposed  acceleration 
increments.  The  answer  in  any  case  depends  on  the  relation  between  wing 
loading  and  mass  parameter.  Figure  39  shows  that  the  substitution  of 
wing  loading  appears  satisfactory  for  transports  but  results  in  too-low 
acceleration  increments  for  personal  airplanes  and  too-high  values  for 
flying  boats. 


4o 


MCA  TN  1976 


Phaae  lag. -  The  data  from  the  general  analysis  were  utilized  to 
obtain  the  lag  in  chords  of  the  acceleration  peak  behind  the  flat-top- 
gust  peak.  Figure  hi  shows  that  the  lag  is  primarily  a  function  of  the 
gradient  distance  for  a  given  mass  parameter  although  some  variations 
due  to  the  tail  area  and  tail  length  are  noted.  No  effect  due  to  a 
variation  in  the  center-of-gravity  position  from  15  to  35  percent  of  the 
chord  is  noted.  The  results  indicate  that  for  a  mass  parameter  of  approxi¬ 
mately  10,  the  lag  in  the  gust  peak  could  be  as  much  as  2  chords  for  an 
8-  to  10-chord  gust.  Computations  indicate  that  this  lag  increases  to 
3  chords  for  a  gradient  distance  of  8  chords  when  the  mass  parameter  is 
about  30. 

In  the  calculations  for  the  triangular  and  sinusoidal  gusts,  the 
lag  of  the  acceleration  peak  behind  the  gust  peak  for  infinite  mass 
parameter,  the  most  adverse  case,  does  not  exceed  2  chords.  The  effect 
of  the  phase  lag  of  the  acceleration  behind  the  gust  can  probably  be 
neglected  since,  for  many  calculations  of  acceleration  increments  and 
load  factors,  a  change  of  1  chord  in  8  has  a  negligible  effect  on  the 
acceleration  ratio. 

Static  margin. -  Figure  34  indicates  that  the  total  wing  load  on  a 
conventional  airplane  decreases  with  increasing  stability  and  increasing 
gradient  distance  because  of  the  resulting  increase  in  pitch.  The 
figure  also  shows  that  the  method  of  obtaining  a  given  static  margin 
is  more  important  than  the  actual  value  of  the  static  margin.  The 
results  for  tailless  airplanes  are  discussed  in  the  following  section. 

Although  no  data  are  available  for  canard  airplanes,  a  comparison 
of  figures  38(b) )  and  32  (both  airplanes  have  the  same  static  margin)  show 
a  marked  effect  on  the  acceleration  increment  due  to  the  radical  change 
in  configuration.  From  the  results  obtained,  it  is  concluded  that  the 
effect  of  airplane  pitch  on  the  wing  load  cannot  be  specified  by  static 
margin  alone. 

Center-of-gravity  position.-  Inspection  of  figure  34  indicates  that, 
for  a  given  configuration,  the  alleviation  due  to  pitch  varies  almost 
directly  with  the  center-of-gravity  position-  The  rate  of  change  of 
pitch  effect  with  center-of-gravity  position  is  roughly  independent  of 
the  configuration  for  the  conditions  analyzed  and  depends  mainly  on  the 
gust-gradient  distance.  For  a  gust  with  a  gradient  distance  of  8  chords, 
figure  34  indicates  that  a  1-percent  change  in  center-of-gravity  posi¬ 
tion  gives  a  1-percent  change  in  load  factor.  The  change  in  load  is 
doubled  for  the  gust  with  a  gradient  distance  of  16  chords  and  about • 
halved  for  a  sharp-edge  gust. 

The  results  shown  in  figure  34,  which  are  for  a  single  value  of  the 
mass  parameter,  and  those  presented  in  figure  36  indicate  that  the  effect 
of  pitch  on  the  gust  load  factor  will  increase  as  the  mass  parameter  is 
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increased.  Figure  36  shows  that  for  a  mass  parameter  of  10,  a  1 -percent 
change  in  center-of -gravity  position  represents  about  a  1-percent  change 
in  load  factor,  whereas  tripling  the  mass  parameter  to  a  value  of  about 
30  increases  the  rate  of  change  in  the  total  wing-load  increment  2  percent 
for  each  1-percent  shift  in  center-of -gravity  position. 

The  tailless  airplane  for  which  calculations  were  made  had  a  mass 

parameter  of  about  26,  and  a  Itt -percent  increase  in  load  for  a  1-percent 

rearward  movement  of  the  center  of  gravity  for  8  chords  was  indicated 
(fig.  38(a)),  which  is  in  fair  agreement  with  the  curves  shown  in  fig¬ 
ure  34  for  the  conventional  airplane.  Although  the  analysis  is  crude, 
the  results  indicate  that,  for  both  the  conventional  and  tailless  air¬ 
planes,  the  effect  of  the  center-of -gravity  position  is  about  the  same. 
Consideration  of  figure  34  indicates  that,  although  the  center  of  gravity 
is  significant  in  determining  the  changes  in  load  on  the  airplane,  the 
center-of -gravity  position  is  not  significant  in  setting  the  absolute 
level  of  the  total  load  for  any  particular  airplane. 

Tail  volume.  -  The  two  middle  lines  of  figure  34  for  approximately 
equal  tail  volume  indicate  that  this  quantity  by  itself  is  not  of  prime 
importance.  Inspection  of  the  results  for  a  gradient  distance  of 
8  chords  shows  that  for  equal  tail  volumes  but  different  tail  lengths  a 
variation  of  10  percent  in  the  effect  of  pitch  on  the  applied  wing  load 
increment  can  be  obtained  for  a  given  static  margin.  For  similar  con¬ 
ditions  and  a  gust  with  a  gradient  distance  of  16  chords,  the  change  in 
load  increment  is  about  20  percent.  A  change  in  tail  volume  obtained 
by  increasing  the  tail  length  and  reducing  the  tail  area  results  in 
higher  loads  on  the  airplane. 

The  results  shown  in  figure  34  for  the  same  static  margin  indicate 
that  the  change  in  tail  area  does  not  yield  significant  variations  in 
load,  whereas  the  change  in  tail  length  results  in  a  significant  modi¬ 
fication  of  the  amount  of  wing  load  imposed  on  the  airplane.  The  data 
thus  indicate  that  for  the  maximum  alleviation  of  load  by  means  of  pitch 
of  the  airplane,  a  short  tail  length  and  large  tail  area  appear  to  be 
most  beneficial,  although  the  combination  will  probably  lead  to  an 
uncomfortable  ride  due  to  the  violent  pitching  motion  of  the  airplane- 
The  use  of  a  small  area  with  a  short  tail  length  will  yield  results  not 
too  different,  although  the  reduction  of  damping  in  pitch  may  emphasize 
the  pitching  oscillation  to  the  detriment  of  passenger  comfort  and 
control  of  the  airplane  in  continuous  rough  air. 

On  the  basis  of  these  analytical  results,  it  is  concluded  that  the 
tail  volume  is  a  secondary  factor  in  establishing  the  level  of  wing  load. 
The  results  indicate  that  the  larger  the  tail  length  for  a  given  static 
margin  the  greater  the  load  imposed  on  the  wing.  For  a  gradient  distance 
of  8  chords  the  effect  on  the  total  wing  load  due  to  the  various  elements 
of  configuration  can  vary  from  10  to  -20  percent. 
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Piloting  and  continuous  rough  air.-  Previous  analytical  and  experi¬ 
mental  studies  were  made  for  a  single  gust  isolated  in  space  and  with 
elevator  fixed,  whereas  airplane  flights  are  in  continuous  rough  air  with 
a  pilot  modifying  the  reactions  of  the  airplane.  Although  there  are  few 
available  data  covering  this  condition,  on  the  basis  of  statistical 
analysis,  significant  differences  have  been  found  to  exist  between  pilots 
and  between  airplanes  of  the  Bame  type.  The  effect  of  different  pilots 
was  as  much  as  20  percent  and  differences  between  airplanes  were  5  to 
10  percent.  It  appears  from  these  results  that  the  pilot  is  more  import¬ 
ant  than  the  variation  between  airplanes  of  the  same  type  in  determining 
the  loads  applied  to  the  airplane-  The  variation  due  to  both  factors  is 
estimated  to  average  15  percent. 

Little  or  no  correlation  between  the  maximum  angular  motions  of  an 
airplane  and  the  maximum  loads  in  the  same  traverse  has  been  found.  It 
was  found  in  another  flight  test  that  one  pilot  moved  the  controls  five 
times  as  often  as  another  pilot  and  no  significant  difference  in  the  leads 
imposed  was  found.  These  results  led  to  the  conclusion  that  the  effect 
of  the  pilot  on  the  imposed  load  is  due  to  his  past  actions  and  is  not 
directly  related  to  the  effect  of  any  single  gust. 

At  the  present  time  the  question  of  the  effect  of  piloting  and  air¬ 
plane  motion  in  continuous  rough  air  on  the  prediction  of  gust  load 
factors  is  only  being  approached.  The  redeeming  feature  at  present  is 
the  fact  that  the  evaluation  of  flight  load  experience  under  actual 
operating  conditions  will  include  piloting  effects  for  the  airplanes 
tested.  The  question  still  to  be  answered  is  whether  data  on  piloting 
effects  taken  in  the  past  on  older  types  of  airplanes  are  valid  for 
modern  high-speed  airplanes  flown  by  pilots  with  different  training  and 
with  airplanes  designed  for  different  requirements  of  flying  qualities 
and  stability. 

TJnsvmmetrical  gusts  and  airplane  response. -  Little  is  known  concern¬ 
ing  the  response  of  airplanes  under  the  action  of  a  gust  which  strikes 
only  one  wing.  The  only  material  available  is  empirical  in  character 
and  has  been  previously  presented  in  this  paper  in  connection  with  the 
studies  of  gust  structure.  In  that  material  it  was  indicated  that  cal¬ 
culation  by  a  simplified  method,  neglecting  the  alleviation  effects  of 
motion,  would  have  to  be  adequate  until  more  information  concerning  the 
significant  parameters  and  the  motions  of  the  airplane  is  obtained. 

Horizontal  tail  loads.-  The  results  given  in  figure  37  indicate  that 
the  static  margin  and  the  tail  volume  are  not  of  primary  importance  in 
the  determination  of  tail  load.  The  data  indicate  a  reasonably  ordered 
variation  of  tail  load  with  center-of-gravity  position.  The  magnitude 
of  the  tail  load  varies  rapidly  with  gradient  distance,  the  variation 
ranging  fham  about  50  percent  for  a  gradient  distance  of  0  chords 
to  0  percent  for  a  gradient  distance  of  16  chords.  Inspection  of 
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t-ables  Xll(a),  Xll(b) ,  and  XII (c)  indicates  that  the  decrease  is  due 
mainly  to  the  increasing  importance  of  the  effect  of  pitch-angle  incre¬ 
ment  as  the  gradient  distance  is  increased  and  due  to  the  fact  that  the 
total  tail  load  is  composed  of  one  large  term  due  to  the  gust,  to  which 
is  added  many  small  contributions  from  various  sources  which  may  be  of 
either  sign.  The  net  result  is  that  the  tail  load  is,  in  many  cases, 
small  and  equal  to  seme  of  the  small  components.  The  calculation  of  the 
tail  load  thus  requires  an  extremely  high  degree  of  accuracy  in  order  to 
obtain  an  accurate  answer. 

As  previously  noted,  the  material  presented  in  figure  37  and  table  XII 
shows  the  tail  load  at  the  time  of  maximum  airplane  acceleration.  It 
is  thought  that  the  maximum  tail  load  may  result  from  a  different 
sized  gust  than  that  which  is  critical  for  the  wing  or  it  may  occur  at  a 
different  time  that  at  maximum  acceleration.  In  view  of  the  many  uncer¬ 
tainties  in  the  calculation  of  tail  loads  due  to  gusts,  it  has  been 
suggested  from  time  to  time  that  the  tail  load  be  calculated  as  the 

increase  in  lift  due  to  the  gust  alone  multiplied  by  the  factor  1  -  4£.- 

da 

A  slightly  different  concept  would  be  to  assume  the  gust  velocity  equal 
to  about  half  that  on  the  wing. 

What  few  experimental  data  from  flight  tests  are  available  appear 
to  be  in  essential  agreement  with  the  preceding  suggestion.  Table  XV" 
indicates  that  the  ratio  of  the  effective  gust  velocity  on  the  tail  as 
determined  from  tail -load  measurements  to  the  effective  gust  velocity 
for  the  wing  is  equal  to  0-53  and  0.64  for  the  XB-15  and  0-2H  airplanes, 
respectively.  The  results  compare  favorably  with  an  estimated  value 

for  1  -  ^  of  0-5- 
da 

In  the  case  of  the  canard  airplane,  inspection  of  the  data  in 
table  XIV  indicates  that  for  a  gust  with  a  gradient  distance  of  8  chords 
the  effects  of  pitch  and  vertical  motion  cancel  so  that  the  total 
stabilizer  load  is  approximately  equal  to  the  load  that  would  be  imposed 
on  the  tail  surface  due  to  the  gust  alone.  For  the  longer  gradient 
distance  of  17*5  chords,  the  cancellation  is  not  complete,  and  the  total 
tail  load  amounts  to  about  two-thirds  of  the  tail  load  resulting  from 
the  action  of  the  gust.  For  the  practical  case,  calculations  of  the 
horizontal  tail  load  for  canard  airplanes  might  be  based  on  the  lift  due 
to  the  change  in  angle  of  attack  resulting  from  the  gust  if  both  the 
vertical  motion  and  pitching  motion  of  the  airplane  are  neglected. 

In  summary,  the  available  results  indicate  that  it  is  not  possible 
by  means  of  detailed  calculations  to  obtain  an  accurate  estimate  of  the 
horizontal  tail  load  due  to  the  action  of  the  gust.  The  tail  load  on 
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the  horizontal  surfaces  can  "be  test  estimated  as  the  tail  load  due  to  the 

gust  alone  multiplied  by  the  factor  1  -  For  the  canard  airplane, 

dot. 

the  downwash  factor  would  be  assumed  to  be  zero. 


Vertical  tail  loads. -  The  suggestion  has  been  made  that  alleviating 
effects  be  neglected  in  the  calculation  of  vertical  tall  loads  due  to 
the  action  of  side  gusts.  On  this  basis,  the  effective  gust  velocity 
for  the  vertical  tail  would  be  related  to  the  vertical  velocity  for  the 
wing  by  the  relation 
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The  experimental  data  in  table  XV  show  excellent  agreement  between  the 
ratio  of  effective  gust  velocities  and  the  reciprocal  of  the  accelera¬ 
tion  ratio.  The  relation  given  is  equivalent  to  the  calculation  of  the 
vertical  tail  load  for  a  true  gust  velocity  with  no  alleviation  due  to 
unsteady-lift  effects  or  airplane  motion-  The  suggestion  appears  satis¬ 
factory  and  it  would  presumably  also  apply  to  the  vertical  tail  on  the 
canard  airplane. 

Steady  lift  in  contrast  to  unsteady  lift. -  Consideration  of  the 
equations  of  motion  for  steady  and  unsteady  flow  indicates  that  the  use 
of  steady-flow  lift  functions  for  generalized  studies  of  gust-load 
problems  is  not  warranted.  The  ratio  of  the  "steady -flow"  acceleration 
to  the  "unsteady -flow"  acceleration  in  the  same  gust  is 
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where  C  is  the  unsteady  lift  coefficient  for  a  gradient  gust.  For 
infinite  mass  parameter,  the  ratio  of  accelerations  varies  from  infinity 
for  zero  gradient  distance  (C  =  0  at  EE  =  0)  to  1.0  for  infinite 
gradient  distance  (steady  flow,  C  =  1.0).  The  variation  with  mass 
parameter  would,  not  he  so  drastic  because  it  affects  only  the  allevia-' 
tion  term,  but,  since  the  integral  in  the  denominator  contains  an 
unsteady-lift  function,  considerable  deviation  would  be  expected  at  small 
values  of  mass  parameter.  Therefore,  the  trends  indicated  by  steady- 
flow  calculations  are  not  necessarily  expected  to  be  those  for  unsteady- 
lift  calculations. 

The  calculated  results  shown  in  figure  38(b) ,  in  which  the  unsteady- 
lift  functions  for  A  =6  and  «  were  utilized,  indicate  the  importance 
of  unsteady -lift- curve  shape  at  large  gust -gradient  distances.  As  indi¬ 
cated  in  the  figure,  for  gradient  distance  of  about  18  chords,  the 
substitution  of  unsteady-lift  functions  of  aspect  ratio  6  for  those  of 
aspect  ratio  00  changed  the  calculated  acceleration  increment  by  50  per¬ 
cent.  It  is  apparent  that  in  the  case  in  which  the  pitching  motion  of 
the  airplane  is  significant  and  the  detailed  response  of  the  airplane  is 
to  be  calculated,  minor  changes  in  the  unsteady-lift  functions  themselves 
can  lead  to  serious  discrepancies.  The  neglect  of  unsteady -lift  relations 
should  involve  more  radical  deviations  than  the  substitution  of  one 
unsteady-lift  function  for  another  and  does  not  appear  to  be  valid. 

Gust  shape . -  Figure  40  shows  that  the  sinusoidal  gust  can  be  sub¬ 
stituted  for  the  triangular  gust  with  a  small  relative  error.  The  curves 
indicate  that  the  relative  error  is  about  2  percent  for  a  range  of  mass 
parameter  from  5  to  80,  with  the  absolute  error  varying  from  5  to 
about  8  percent  high.  The  results  indicate,  therefore,  that  the  trends 
predicted  for  the  one  gust  shape  should  apply  equally  well  to  the  other 
gust  shape. 

Calculated  and  experimental  results. -  Inspection  of  figure  28 
indicates  that,  for  the  case  in  which  the  pitch  effects  are  negligible, 
the  scatter  in  computing  the  acceleration  ratio  appears  to  be  random 
and  the  experimental  results  fall  on  each  side  of  the  calculated  curve. 
Some  of  the  scatter  noted  is  caused  by  pitching  of  the  model  and  some 
is  undoubtedly  caused  by  difficulties  in  making  the  experimental  measure¬ 
ments.  On  the  whole,  the  calculated  results  appear  to  agree  in  a  satis¬ 
factory  manner  with  experiment.  Similar  agreement  is  indicated  in 
figures  35>  38(a),  38(b),  and  32  for  H  =  0  where  pitch  is  considered 
negligible.  It  Is  believed,  therefore,  on  the  basis  of  these  results, 
that  the  unsteady-lift  functions,  the  slope  of  the  lift  curve,  and 
similar  items  involved  in  the  calculations  are  known  with  an  accuracy 
comparable  to  that  for  other  elements  of  gust-load  calculations. 
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For  gradient  distances  other  than  zero  or  for  those  cases  in  which 
the  pitch  is  a  factor,  the  discrepancies  "between  calculation  and  experi¬ 
ment  "become  larger  as  the  gradient  distance  increases  until,  as  indicated 
"by  figure  35,  the  errors  are  serious  for  gradient  distances  of  16  to 
chords.  In  the  case  of  "both  the  conventional  airplanes  (fig.  35)  and 
the  canard  airplane  (fig-  38(b) ),  the  calculations  for  the  longer  gust- 
gradient  distance  are  unconservative  and  lead  to  smaller  predicted  load 
factors.  In  the  case  of  the  tailless  airplane  (fig.  38(a))  the  calculated 
values  are  in  good  agreement.  Consideration  of  the  data  presented  in 
figure  32,  in  which  the  effect  of  pitch  has  "been  neglected  in  the  cal¬ 
culations,  indicates  that  the  neglect  of  pitch  for  the  longer  gust- 
gradient  distances  will  "be  conservative.  The  good  agreement  for  the 
tailless  airplane  (fig.  38(a))  at  all  gradient  distances  compared  to  the 
variation  in  agreement  "between  calculation  and  experiment  for  the  other 
cases  indicates  that  the  introduction  of  the  tail  surface  and  possibly 
the  fuselage  has  much  to  do  with  the  adequacy  of  the  detailed  calcula¬ 
tions  of  airplane  response  to  a  gust.  The  errors  indicated  for  the 
longer  gust-gradient  distances  as  compared  to  the  sharper  gust  indicate 
that  the  pitch  effect  is  probably  the  least  accurately  predicted  quan¬ 
tity  and,  therefore,  leads  to  the  errors  noted. 

From  the  results  obtained  it  appears  that  the  wing  load  due  to  a 
gust  can  be  calculated  within  10  percent  for  gradient  distances  up  to 
10  chords.  The  accuracy  in  any  individual  case  or  for  large  gradient 
distances  cannot  be  estimated  because  of  pitch  and  secondary  effects 
that  may  be  significant  but  are  not  recognized  until  test  results  are 
available.  Minor  changes  in  variables  such  as  downwash  can  lead  to 
serious  variation  of  numerical  results,  and  the  estimation  of  loads  by 
detailed  or  extensive  calculations  is  not  of  sufficient  accuracy  to 
warrant  much  confidence  in  the  results. 


ELASTIC -AIRPLANE  REACTIONS 


Since  the  loads  imposed  when  airplanes  encounter  gusts  are  applied 
suddenly,  the  dynamic  response  of  the  airplane  structure  has  been  of 
concern  since  the  initiation  of  gust -load  studies.  A  number  of  studies 
(references  2  and  35  to  40)  have  been  made  at  various  times  to  evaluate 
the  importance  of  dynamic  response  and  of  the  various  parameters  involved. 
Icl  1939,  projected  airplane  designs  indicated  that  dynamic  response  might 
be  of  concern  and  an  analytical  and  experimental  study  was  undertaken. 

The  results  of  this  investigation  are  presented  in  reference  38. 

In  this  section  the  results  obtained  in  reference  38  are  summarized 
as  well  as  some  results  from  unpublished  studies.  The  significance  of 
dynamic  response  and  the  importance  of  the  various  parameters  in  producing 
dynamic  responses  in  airplane  structures  are  of  chief  concern. 
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Methods 

The  analytical  method,  developed  in  reference  38  was  to  reduce  the 
airplane  structure  to  an  "equivalent"  biplane  whose  upper  wing  had  the 
same  motion  as  the  original  wing  tip  and,  hy  using  an  effective  damping 
factor  .instead  of  unsteady -lift  functions,  to  obtain  two  simultaneous 
linear  differential  equations-  The  equivalent  biplane,  in  which  the 
upper  wing  is  connected  by  springs  to  the  lower  wing-fuselage  combina¬ 
tion,  is  adjusted  so  that  the  components  have  the  same  motions  as  the 
wing  tip  and  fuselage  of  the  airplane  under  study.  The  equivalent 
system  must  include  the  proper  distribution  of  aerodynamic  forces  as 
well  as  inertia  and  elastic  forces.  The  deflection  of  the  upper  wing 
of  the  biplane  is  taken  as  a  measure  of  dynamic  stress.  The  dynamic- 
stress  ratio  is  defined  as  the  ratio  of  dynamic  to  static  wing-tip 
deflections-  The  static  deflection  is  computed  in  the  same  manner  as 
for  normal  design  by  including  inertia  effects  but  neglecting  aerodynamic 
damping  due  to  vibration  of  the  upper  wing. 

The  spring  constant,  equivalent  masses,  and  aerodynamic  damping  for 
the  equivalent  system  are  calculated  so  that  the  static  and  vibration 
characteristics  of  the  original  wing  are  represented.  In  general,  the 
following  conditions  are  to  be  satisfied: 

1.  The  total  mass  of  and  the  total  load  on  the  equivalent  biplane 
should  be  identical  with  those  of  the  original  airplane. 

2.  The  upper  wing  should  deflect  under  the  equivalent  static  load 
the  same  amount  as  the  original  wing  tip  under  its  corresponding  static¬ 
load  distribution. 

3.  The  natural  frequency  of  the  equivalent  system  should  be  the 
same  as  that  of  the  original  wing. 

h.  The  kinetic  energy  of  vibration  of  the  upper  wing  should  closely 
approximate  that  of  the  original  wing  for  the  same  tip  amplitude  of 
vibration. 

5-  The  damping  coefficient  of  the  upper  wing  should  represent,  at 
least  up  to  peak  load,  the  damping  of  the  motion  of  the  original  wing. 

The  shape  of  the  forcing  function  used  in  the  calculations  of 
reference  38  was  obtained  from  accelerometer  records  of  gust-tunnel  tests 
of  different  models  by  subtracting  the  computed  acceleration  increment 
due  to  vertical  motion.  This  procedure  was  followed  since  the  prediction 
of  airplane  reactions  as  affected  by  stability  and  other  factors  is  of 
doubtful  accuracy,  as  previously  noted.  For  representative  gust  sizes 
(H  =  0  to  H  =  20  chords)  it  was  found  in  reference  38  that  a  curve  of 

the  form  Ate-"^  was  a  good  representation.  The  forcing  function  for 
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the  Complete  airplane  was  then  divided  so  that  the  upper  wing  of  the 
equivalent  "biplane  would  have  the  same  static  deflection  as  the  original 
wing  tip. 

In  addition  to  and  as  a  check  on  the  method  of  calculation,  tests 
were  made  in  the  Langley  gust  tunnel  with  a  model  having  two  wings  of 
different  frequencies.  The  wing  deflection  and  fuselage  acceleration 
were  the  primary  quantities  measured.  The  model  used  is  shown  diagram- 
matically  in  figure  42  and  its  characteristics  together  with  the  test 
conditions  are  given  in  table  XVI. 

For  each  of  the  two  wing  frequencies  tests  were  made  at  one  forward 
speed  and  three  gust  sizes.  Time  histories  of  pitch,  acceleration  incre¬ 
ment  of  the  fuselage,  and  wing  deflection  were  obtained.  Some  results 
of  the  tests  are  shown  in  figure  43,  in  which  the  maximum  wing-tip  deflec¬ 
tion  per  unit  acceleration  increment  is  plotted  for  one  wing  as  a  func¬ 
tion  of  gradient  distance. 

The  analytical  procedure  described  in  reference  38  was  utilized  to 
compute  the  response  of  the  system  when  damping  of  the  motion  is  included 
and  when  aerodynamic  damping  is  neglected.  The  results  have  been  included 
in  figure  43  for  comparison  with  experiment. 


Analytical  Study 

In  order  to  determine  the  effect  of  the  various  parameters  on 
dynamic  response,  a  series  of  calculations  was  made  at  three  gust- 
gradient  distances  for  four  airplanes.  The  characteristics  of  the  air¬ 
planes,  which  were  designated  as  A,  B,  C,  and  D,  are  given  in  table  XVII. 
The  parameters  studied  were  the  wing  stiffness,  airplane  weight,  and 
forward  speed,  with  all  damping  neglected.  The  gust  velocity  was  not 
varied  since  ratios  of  dynamic  to  static  conditions  were  utilized  and  it 
was  not  a  factor  for  linear  equations.  The  conditions  used  in  the  cal¬ 
culations  are  given  in  table  XVII* 

The  calculated  results  for  each  airplane  of  the  maximum  deflection 
ratio,  wing-tip  acceleration  ratio,  and  fuselage  acceleration  ratio  are 
shown  in  figure  44. 


Discussion 

Inspection  of  the  results  of  figure  43  indicates,  in  general,  that 
dynamic  response  is  overestimated  when  damping  is  neglected;  whereas  the 
inclusion  of  damping  leads  to  fair  agreement  with  experimental  results. 

The  results  obtained  were  not  sufficiently  accurate  to  provide  an  absolute 
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check  of  the  calculation.  The  results  do  indicate,  however,  that  the 
calculations  give  a  fair  prediction  of  the  ratio  of  maximum  wing  deflec¬ 
tion  to  maximum  fuselage  acceleration. 

Time  histories  given  in  reference  38  "but  not  shown  herein  indicate 
that  the  method  of  calculation  did  not  predict  correctly  the  amplitude 
of  the  wing  oscillations  after  passing  the  peak  load.  The  result  is 
due  in  part  to  the  fact  that  the  damping  coefficient  was  adjusted  to 
give  a  good  estimate  for  the  first  wing  motion  and  would  he  overestimated 
for  the  subsequent  vibratory  motion.  Putnam  (reference  39)  attempted 
to  resolve  this  difficulty  when  he  extended  the  method  of  reference  38. 
Another  possible  factor  is  that  the  simplifications  used  are  still  too 
drastic.  Other  results  (reference  38)  indicate  that  the  use  of  a  constant 
damping  factor  is  reasonable  if  its  magnitude  is  adjusted  for  the  average 
effect  of  unsteady  lift. 


Inspection  of  figure  44  indicates  that  the  dynamic -stress 


ratios 


'Vax/5 


and  wing-tip  acceleration  ratios 


increase  as  the  gradient  distance  decreases-  Recent  calculations  for 
more  modern  aircraft  bear  out  this  result  and  indicate  that  elastic 


response  is  becoming  increasingly  important.  The  fuselage  acceleration 

ratio  D2&.P  /An~.  did  not  appear  to  differ  greatly  from  1-0  and 

max/  max 


does  not  appear  to  be  seriously  affected  by  gust  size. 


The  calculations  shown  in  figure  45  for  the  effect  of  wing  stiffness 
on  the  dynamic -stress  ratios  indicate  that  at  the  "critical"  gradient 
distance  of  10  chords,  the  high-frequency  wing  shows  a  dynamic -stress 
ratio  of  14  percent  below  the  low-frequency  wing.  Further  analysis  is 
needed  before  a  definite  conclusion  can  be  reached  that  a  reduction  of 
wing  frequency  by  changing  wing  stiffness  tends  to  increase  the  dynamic- 
stress  ratio  at  any  gradient  distance. 


The  results  given  in  figure  46  illustrate  the  variations  of  the 
three  ratios  due  to  increasing  the  forward  velocity  of  model  C  from 
200  to  400  miles  per  hour.  The  increase  in  velocity  together  with  the 
corresponding  increase  in  the  rate  of  application  of  the  gust  load  would 
appear  to  result  in  an  increase  in  the  dynamic -stress  ratio.  Figure  46, 
however,  shows  that  the  dynamic -stress  ratio  does  not  vary  much  as  the 
speed  increases  and  this  lack  of  variation  is  thought  to  be  caused  partly 
by  the  increase  in  aerodynamic  damping  with  speed-  Although  the  results 
for  the  fuselage  acceleration  ratio  show  little  effect  of  speed,  the 
wing-tip  acceleration  ratio  increases  from  about  1.8  to  2.5  as  the  speed 
is  doubled. 


The  calculations  for  model  D  were  made  to  show  the  effect  on  the 
dynamic-stress  ratio  of  a  change  in  flight  condition  from  normal  gross 
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weight  to  overload  gross  weight.  Table  XVII  shows  that  the  forward 
velocity  is  different  in  the  two  cases,  "but  consideration  of  the  fore¬ 
going  discussion  may  Justify  the  assumption  that  speed  has  a  negligible 
effect  on  the  findings.  The  results  are  given  in  figure  47  with  dynamic- 
stress  ratio  plotted  as  a  function  of  the  gradient  distance  of  the  gust. 

A  reduction  in  wing  frequency  brought  about  by  the  addition  of  mass  is 
shown  to  result  in  an  increase  in  the  dynamic -stress  ratio. 

Available  information  on  repeated  gusts  such  as  given  in  table  IV 
has  been  used  (reference  38)  to  determine  the  dynamic  effects  on  wing 
deflection  in  repeated  gusts.  The  following  table  (from  reference  38) 
indicates  the  dynamic -deflection  ratios  for  two  gusts  of  given  intensity 
spaced  25  chords  apart.  Also  given  in  the  table  are  the  deflection 
ratios  for  the  single  gusts  of  greater  intensity  which  have  the  same 
probability  of  occurrence  as  the  sequence  of  two  gusts  used. 


Model 

Condition 

II—  1  11 

Design 

repeated 

gusts 

(5<w/5st) 

Design 
single 
gusts  > 

(BW5st) 

C 

1 

0.96 

1.07 

D 

1 

1.08 

.92 

D 

2 

1.10 

1.09 

The  variation  in  the  results  indicates  that  the  dynamic-stress  ratios 
for  a  repeated  gust  should  be  investigated.  The  results  also  indicate 
that  these  dynamic -stress  ratios  are  probably  not  much  greater  than 
those  determined  for  a  single  gust  likely  to  be  encountered  in  flight. 


Concluding  Remarks  Concerning  Elastic -Airplane  Reactions 

Study  of  available  results  and  unpublished  data  indicates  that 
dynamic  response  is  becoming  of  greater  importance  with  modern  advances 
in  airplane  design.  As  mass  tends  to  be  distributed  along  the  wing, 
more  exact  solutions  are  needed  and,  until  the  various  factors  such  as 
unsteady-lift  functions,  inertia  terms,  elastic  constants,  and  the  span- 
wise  gust  distribution  are  known  more  exactly,  such  calculations  should 
■be  utilized  on  a  relative  basis.  Although  new  and  more  elaborate  methods 
of  calculation  are  available,  serious  problems  still  remain  as  to  the 
basic  aerodynamic  and  elastic  coefficients. 
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OPERATING  STATISTICS 


■  A  knowledge  of  gust  structure  and  methods  for  computing  airplane 
reactions  is  insufficient  for  gust -load  calculations  without  knowing' the' 
conditions  for  which  load  calculations  are  required.  Thus,  in  order  to 
solve  the  practical  problem  of  load  prediction,  the  gusts  that  airplanes 
encounter  and  the  operating  conditions  (speed,  weight,  altitude,  and, 
perhaps,  center-of -gravity  position)  that  are  likely  to  exist  at  the 
time  the  gust  is  encountered  must  he  known.  The  fact  that  the  source  of 
load,  the  atmospheric  gust,  is  of  a  fairly  random  character  and  that  the 
operating  conditions  vary  widely  make  it  impractical  to  predict  the 
exact  loads  and  the  associated  conditions  that  are  experienced  by  a 
given  airplane  during  its  lifetime. 


METHOD 


The  general  method  of  obtaining  the  desired  information  has  been 
to  install  instruments  in  aircraft  and  to  record  their  experiences  and 
operating  conditions  without  interfering  with  routine  practices.  Three 
procedures  have  been  used:  the  installation  of  simple  instruments  such 
as  the  NACA  V-G  recorder  giving  broad  coverage  as  to  route,  airplane 
type,  and  time  but  no  detailed  information,  the  installation  of  special 
instrumentation  for  a  very  limited  period  of  time  to  obtain  information 
on  a  particular  quantity  in  detail,  and,  finally,  the  installation  of 
fairly  elaborate  instrumentation  controlled  by  an  observer  to  obtain 
as  much  detailed  data  as  possible  on  several  quantities  during  one  or 
more  flights. 

The  NACA  V-G  recorder  (reference  4)  records  on  a  smoked  glass  plate 
the  normal  acceleration  as  a  function  of  the  airspeed  at  which  it  was 
imposed.  In  use,  a  record  plate  is  left  in  the  instrument  during 
operations  and  the  resulting  clear  area  on  the  glass  represents  an 
envelope  of  accelerations  and  airspeeds  experienced  by  the  airplane. 

Standard  NACA  photographically  recording  instruments  or  motion 
pictures  of  indicating  instruments  are  used  for  more  detailed  measurements. 


SCOPE  OF  DATA 


The  characteristics  and  pertinent  dimensions  of  transport  airplanes 
for  which  gust  statistics  have  been  collected  are  given  in  table  XVIII. 

The  listed  values  of  the  limit  load-factor  increment  due  to-  gusts,  the  high 
speed  in  level  flight,  and  the  placard  do-not -exceed  speed  were  obtained 
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when  possible  from  official  sources.  For  the  older  airplanes  where  the 
design  conditions  differed  from  modern  requirements,  the  pertinent  speeds 
and  load  factors  were  recomputed  according  to  reference  1  in  order  to 
place  all  the  characteristics  on  a  comparable  basis. 

Table  XIX  is  a  summary  of  the  conditions  investigated  and  includes 
the  routes,  periods  covered,  and  the  amount  of  data  obtained.  The  opera¬ 
tions  have  been  subdivided  into  prewar,  wartime,  and  postwar  periods. 

For  the  prewar  and  wartime  periods,  information  was  obtained  with  the  V-G 
recorder,  but  for  the  postwar  period  sufficient  V-G  records  are  not 
available.  Special  studies  of  speed  and  Mach  number  variations  are 
available  and  some  of  the  results  have  been  included.  The  material 
covered  in  table  XIX  does  not  include  all  the  data  obtained  since  scat¬ 
tered  V-G  data  of  insufficient  scope  for  analysis  have  been  disregarded. 


STATISTICAL  METHODS 


The  random  character  of  gusts  and  the  lack  of  control  over  flight 
conditions  have  resulted  in  the  use  of  statistical  methods  of  analysis 
to  smooth  data,  to  eliminate  improper  weighting,  to  extrapolate  results, 
and  to  place  data  From  different  sources  on  a  comparable  basis.  The 
basic  data  in  most  cases  are  a  count  of  the  values  of  a  quantity  according 
to  magnitude.  Since  the  application  of  statistical  methods  to  gust  loads 
is  fairly  recent,  the  mass  of  data  collected  earlier  offers  some  diffi¬ 
culties  because  the  requirements  of  statistical  analysis  were  not 
considered  in  their  collection. 

Pearson  type  III  probability  distribution  curves  (see  references  hi 
and  42)  have  been  utilized  in  the  analysis  on  the  assumption  that  they 
adequately  represent  the  data  on  operating  statistics.  The  curve  is 
determined  by  three  parameters:  the  mean  value,  the  standard  deviation, 
and  the  coefficient  of  skewness.  In  general,  the  goodness  of  fit  has 
been  based  on  engineering  Judgment  rather  than  any  test  procedure- 

An  important  problem  in  the  study  of  the  frequency  of  exceeding  the 
larger  values  of  load  or  speed  -under  operating  conditions  is  the  deter¬ 
mination  of  whether  observed  or  estimated  differences  in  frequencies 
between  samples  are  real  or  represent  limitations  of  the  data.  No 
satisfactory  test  for  this  purpose  has  been  found,  but  a  5;1  ratio  of 
frequencies  has  been  used  in  reference  43  as  an  engineering  measure.  If 
the  differences  are  less  than  this  criterion,  real  differences  may  exist 
but  are  considered  too  small  to  be  determined  from  the  available  data. 

A  serious  limitation  in  the  analysis  of  operating  statistics  exists 
since  the  results  are  used  to  predict  future  operating  experience.  If 
no  extraneous  factor  enters,  such  as  changes  In  design  rules,  operating 
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regulations,  or  the  introduction  of  never  airplanes,  then  the  prediction 
may  he  satisfactory.  Such  changes  as,  for  example,  the  shift  from  pre- 
var  to  postwar  operations  may  lead  to  some  uncertainty  as  to  the  appli¬ 
cability  of  a  prediction.  Predictions  assume  that  the  trend  in  the  data 
has  not  been  affected  by  new  factors.  (No  method  yet  is  available  for 
evaluating  or  detecting  the  effect  of  changes  in  the  data  now  used. ) 

Other  questions  that  occur  in  statistical  analyses  arise  in  connec¬ 
tion  with  the  existence  or  inclusion  of  data  from  two  different  regions 
or  populations  as  part  of  the  same  sample,  the  presence  of  absolute  and 
physical  boundaries  that  cannot  be  exceeded,  the  effect  of  extraneous 
variables,  and  the  effect  of  gradual  transitions  from  one  region  to 
another.  The  presence  of  such  limitations  can  have  serious  results  as 
far  as  the  significance  of  the  results  is  concerned  and  such  effects 
cannot  always  be  evaluated  or  separated.  An  example  of  the  effect  of 
different  regions  might  be  the  combination  of  acceleration  data  involving 
flight  below  maximum  lift  coefficient  and  flight  above  maximum  lift 
coefficient.  The  aerodynamic  characteristics  of  the  airplane  and  its 
behavior  in  the  two  regions  are  entirely  different.  Thus,  a  sample  that 
combines  such  data  does  not  offer  an  adequate  means  of  predicting  future 
expectations  of  a  given  acceleration  except  for  the  conditions  of  the 
original  data,  since  any  derived  curves  are  a  composite  of  two  independent 
sets  and  therefore  the  prediction  can  only  apply  to  the  same  combination. 


RESULTS 


V-G  data.  -  From  the  individual  V-G  records  that  are  summarized  in 
table  XIX,  frequency  distributions  of  the  maximum  acceleration 
increment  An^^,  the  maximum  speed  V  Y,  and  the  speed  VQ  at  which 

maximum  acceleration  increment  was  experienced  have  been  compiled  for 

the  wartime  and  prewar  periods.  The  statistical  parameters:  the  mean, 

the  standard  deviation,  and  the  skewness  are  given  in  table  XX,  which 

also  includes  the  values  of  V0  ,  An  ,  and  V  recorded  during 

max7  max7  max  ^ 

the  period  under  consideration. 

The  parameters  of  the  Pearson  type  III  curves  given  in  table  XX 
have  been  used  to  obtain  the  flight  miles  to  exceed  the  specified  values 
of  a  selected  quantity.  The  transformation  from  probability  to  flight 
miles  was  performed  through  use  of  a  nominal  cruising  speed  and  the 
average  flight  hours  per  V-G  record.  Table  XXI  gives  the  flight  miles 
to  exceed:  (l)  the  limit  load-factor  increment,  (2)  the  placard  speed  of 
the  airplane,  and  (3)  an  acceleration  increment  corresponding  to 
encountering  a  gust  with  an  effective  gust  velocity  of  37* 5K  feet  per 
second  at  the  most  probable  value  of  VQ,  Vp.  The  probable  speed  Vp 
has  been  included  in  table  XXI  as  a  fraction  of  the  high  speed  in 
level  flight  V-^- 
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Time -hi story  data.  -  Time  histories  from  the  special  investigations 
of  the  postwar  period  have  teen  evaluated  to  obtain  the  flight  miles  to 
exceed  the  placard  speed,  the  probable  speed  of  flight  for  each  sample, 
and  the  distribution  of  airspeed.  Because  of  the  limited  duration  of 
the  samples,  each  sample  was  read  to  obtain  the  maximum  speed  during  a 
fixed  interval  of  time  (6  to  10  minutes)  and  the  resulting  frequency 
distributions  were  then  used  to  obtain  the  flight  miles  to  exceed  the 
placard  speed.  The  results  are  included  in  table  XXI  and  in  figure  48. 
Figure  48  presents  the  percent  of  total  flight  time  spent  at  speeds 
equal  to  or  greater  than  any  selected  fraction  of  the  high  speed  V^- 

The  probable  speed  noted  in  table  XXI  for  airplanes  E  and  F  is  the 
probable  speed  of  flight  and  is  not,  as  in  the  case  of  other  airplanes, 
the  probable  speed  at  which  maximum  acceleration  will  occur.  A  simple 
analysis  in  which  the  velocity -acceleration  envelopes  were  synthesized 
from  the  speed-frequency-distribution  data  and  the  gust -frequency  data 
of  reference  9  indicated  that  the  probable  speed  for  maximum  accelera¬ 
tion  was  about  0.037-^  greater  than  the  probable  speed  of  flight. 

The  data  obtained  have  also  been  evaluated  according  to  flight 
condition  to  obtain  curves  of  the  flight  hours  required  in  climb,  cruise, 
and  descent  to  exceed  given  values  of  speed.  Figure  k9  is  typical  of  the 
results  obtained  for  airplane  E. 

I 

Disturbed  motions.-  from  the  data  obtained  in  the  special  investi¬ 
gations  a  statistical  study  has  been  made  of  other  quantities  of  interest 
in  gust-load  studies,  such  as  the  effects  of  disturbances  on  the  relative 
frequency  distribution  and  the  fraction  of  the  time  spent  in  rough  air. 

Some  question  has  arisen  as  to  the  effect  of  the  selected  datum  on 
the  frequency  distributions  obtained  in  continuous  rough  air.  This 
question  is  partly  answered  by  the  tests  reported  in  reference  9  on  one 
of  the  roughest  flights  with  the  XC-35  airplane  where  the  relative  gust- 
frequency  distribution  was  determined  by  using  a  disturbed  datum  and  an  ■ 
arbitrary  1  g  datum-  The  results  of  the  study  are  given  in  figure  50* 

Camera  records  of  the  pilot's  instrument  panel  in  the  XC-35  airplane 
have  been  evaluated  to  obtain  the  relative  frequency  of  occurrence  of  the 
maximum  total  variations  in  the  angular  displacements  of  the  XC-35  air¬ 
plane  during  separate  traverses  through  clouds.  The  results  are 
summarized  in  figure  51  as  the  relative  frequency  of  exceeding  selected 
values  of  yawing,  rolling,  and  the  pitching  displacements  and  the  rate 
of  turn.  It  should  be  noted  in  considering  this  figure  that  the  quantity 
plotted  is  the  sum  of  the  maximum  positive  and  the  maximum  negative 
values  recorded  during  each  cloud  traverse. 

Path  ratio. -  The  available  data  for  prewar  transport  operation 
given  in  reference  9  indicate  that  the  path  ratio  (the  miles  of  rough 
air  divided  by  the  total  miles  flown)  varied  from  0.2h  to  0-006  for 
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different  routes  in  the  U*S.  and  abroad  and  had  an  average  value  of  0.1. 
Since,  on  the  average,  500/c  gusts  with  an  effective  gust  velocity  greater 
than  0. 3  foot  per  second  are  encountered  per  mile  of  rough  air  (refer¬ 
ence  9) ,  the  number  of  gusts  that  will  he  encountered  on  the  average 
during  the  lifetime  of  an  airplane  can  he  estimated.  Wo  path-ratio 
data  are  available  for  postwar  operations.  Attempts  to  obtain  such  data 
from  the  indirect  evidence  of  V-G  records  lead  to  erroneous  results. 


DISCUSSION 


Applied  acceleration  increments. -  For  the  prewar  period,  table  XXI 
shows  a  wide  scatter  in  the  flight  miles  to  exceed  the  limit  accelera¬ 
tion  increment.  The  values  vary  by  a  factor  of  about  300.  A  variation 
by  a  factor  of  42  occurs  for  the  same  type  of  airplane  operated  on 
different  routes  by  different  airlines.  In  four  out  of  the  six  samples, 
the  flight  miles  to  exceed  the  limit  acceleration  increment  would  be 
greater  than  13  million,  or,  if  a  cruising  speed  of  200  miles  per  hour 
is  assumed,  it  might  be  expected  that  the  limit  load  factor  would  be 
exceeded  about  once  on  the  average  in  every  6 0,000  hours  of  flight. 

The  data  for  the  wartime  period  do  not  show  as  much  variation  of 
flight  miles  and  the  variations  do  not  appear  significant  within  that 
group.  In  comparison  with  the  prewar  operations,  however,  the  flight 
miles  to  exceed  limit  acceleration  have  consistently  decreased^  this 
fact  indicates  that  the  pressure  of  the  emergency  on  wartime  commercial 
operations  resulted  in  higher  imposed  loads  on  the  airplane. 

Atmospheric  gustiness. -  The  flight  miles  to  exceed  the  acceleration 
increment  corresponding  to  an  effective  gust  velocity  of  37- 5K  feet  per 
second  at  the  probable  speed  (table  XXl)  indicate  that  the  operational 
experience  on  the  basis  of  a  "gust  intensity"  have  much  less  scatter. 

The  spread  in  flight  miles  was  15:1  for  the  prewar  operations  and 
about  7:1  for  the  wartime  period.  For  prewar  conditions,  when  the 
data  for  airplane  D  on  route  V  are  neglected,  the  scatter  is  within  the 
arbitrary  criterion  and,  in  comparison  with  the  spread  in  flight 
miles  for  limit  acceleration  increment,  indicates  that  the  difference 
in  the  level  of  roughness  encountered  on  various  routes  is  not  of 
engineering  concern.  A  similar  observation  is  indicated  for  the  wartime 
period  although  the  smaller  average  flight  miles  indicates  that  flights 
during  that  period  were  through  more  severe  weather  conditions  than  for 
the  prewar  period. 

In  connection  with  prewar  data  for  airplane  D  on  route  V  of 
table  XXI,  the  low  value  of  flight  miles  shown  may  be  significant  and 
indicates  that  early  transpacific  flights  encountered  more  severe 
weather  than  was  experienced  along  other  routes.  When  the  distances 
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involved  in  transpacific  operations  and  the  lack  of  weather  ships  in  this 
period  are  considered,  more  frequent  accidental  encounters  with  severe 
weather  might  he  expected  than  for  the  transcontinental  routes  or 
operations  in  more  populated  regions. 

On  the  basis  of  the  information  concerning  the  flight  miles  to 
equal  or  exceed  the  acceleration  increment  corresponding  to  37- 5K  feet 
per  second  at  the  probable  speed,  the  level  of  route  roughness  is  con¬ 
cluded  to  be  largely  independent  of  the  route,  airplane,  or  operator. 

It  is  not  possible  at  this  time  to  state  whether  the  maintenance  of  a 
constant  level  of  roughness  is  due  to  dispatching  practices,  meteorolo¬ 
gical  forecasting  abilities,  or  both. 

Frequency  of  encountering  gusts. -  The  available  data  on  path  ratio 
given  in  reference  9  indicate  a  wide  spread  in  values  obtained,  and  at 
this  time  an  average  value  of  0.1  seems  to  be  the  best  estimate  avail¬ 
able.  Since  the  path  ratio  is  the  proportion  of  the  total  miles  flown 
that  are  spent  in  rough  air  and  the  number  of  gusts  per  mile  of  rough 
air  has  been  found  to  be  essentially  constant  at  500/c,  the  path  ratio 
can  be  defined  as  the  actual  number  of  gusts  divided  by  the  total  number 
expected  if  the  total  flight  path  were  rough.  It  should  be  noted  that 
the  path  ratio  defined  by  gust  counts  depends  on  the  threshold  (in  this 
case,  the  threshold  is  0.3  fps).  When  used  with  a  gust  count  of 
500/c  gusts  per  mile  of  rough  air  and  the  gust  frequency  distributions  of 
figure  8,  the  recommended  average  value  of  path  ratio  can  be  used  to 
estimate  the  probability  of  an  airplane  encountering  a  gust  of  any 
specified  intensity. 

Although  no  information  is  available  as  to  the  variation  of  the 
number  of  gusts  or  path  ratio  with  altitude,  terrain,  or  weather,  the 
path  ratio  would  probably  decrease  with  altitude. 

Probable  speed  Vp. -  The  determination  of  the  miles  to  exceed  37- 5K 

feet  per  second  at  the  probable  speed  showed  the  flight  miles  to  be 
relatively  constant,  and,  therefore,  the  wide  variations  in  the  flight 
miles  for  limit  load  factor  might  be  ascribed  to  variations  in  the 
probable  speed.  On  this  basis,  the  results  indicate  that  the  probable 
speed  is  Important  in  determining  flight  miles  to  limit  acceleration 
Increment . 

The  data  on  the  probable-speed  ratio  (table  XXl)  indicate  a 
scatter  in  the  speed  ratios  for  any  given  period.  The  speed  ratio  for 
maximum  acceleration  increased  from  the  prewar  period  to  the  postwar 
period  for  the  same  airplanes  and  routes  from  an  average  value  of  0-75 
to  about  0-86.  It  appears  that  the  speed  ratio  might  be  a  function  of 
the  route,  airline  policy,  and  the  airplane  characteristics,  but  no 
conclusion  can  be  drawn  as  to  the  significant  parameters  that  determine 
the  speed  ratio. 
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Although  no  data  have  "been  analyzed  with  regard  to  the  imposed 
gust  loads  for  the  postwar  period,  on  the  basis  of  the  current  design 
requirements  and  if  the  level  of  roughness  is  assumed  to  remain  the 
same  as  for  the  prewar  period,  the  imposed  gust  acceleration  increment 
for  a  given  number  of  flight  hours  will  he  about  12  percent  higher  on 
the  average  for  the  postwar  period  than  for  the  prewar  period-  Alter¬ 
natively,  the  increased  probable  speed  for  the  postwar  period  would 
result  in  an  appreciable  decrease  in  the  flight  miles  to  exceed  a  given 
acceleration  increment- 

Maximum  speeds- -  Inspection  of  figure  49  indicates  that  the 
probability  of  exceeding  a  given  value  of  airspeed  is  a  function  of  the 
flight  condition.  Figure  k9  indicates  that  the  descent  is  most  import¬ 
ant,  and  for  one  postwar  airplane  the  placard  speed  might  be  exceeded 
on  the  average  about  once  in  200  hours  in  descent.  A  simple  sorting  of 
the  data  with  regard  to  the  smooth  air  and  rough  air  showed  no  apparent 
effect  of  rough  air  on  the  speed  tendencies. 

For  complete  flight  operations,  the  data  given  in  table  XXI  show 
that  the  flight  miles  to  exceed  the  placard  speed  decreased  from  many 
millions  of  miles  during  the  prewar  period  to  about  a  quarter  of  million 
of  miles  during  the  postwar  period.  The  reduction  in  the  flight  miles 
to  exceed  the  placard  speed  for  the  postwar  period  is  of  considerable 
significance. 

Sneed-time  distributions.-  Figure  M3  shows  that  a  relatively  higher 
percentage  of  time  is  spent  above  high  speed  in  level  flight  than  would 
be  expected  from  estimates  made  by  airline  personnel,  which  estimates 
have  indicated  a  value  of  less  than  1  percent.  The  high  percentage  may 
be  due  in  part  to  performance  abilities  of  the  airplanes  in  excess  of 
that  indicated  by  the  regulated  speeds  or  due  to  high-speed  descents  at 
the  end  of  flights.  Until  further  information  is  obtained,  the  results 
in  figure  M3  might  be  assumed  to  be  representative  of  postwar  operations. 

Disturbed  motions. -  Tolefson  (reference  44)  showed  that  the  maximum 
angular  displacements  of  an  airplane  do  not  correlate  with  the  maximum 
intensity  of  gusts  encountered  during  a  given  run.  The  correlation 
coefficient  between  the  angular  motions  and  the  maximum  gust  intensity 
was  about  0.  3-  The  disturbed  motions  may  depend  more  on  the  sequence  of 
the  gusts  encountered  than  on  the  intensity  cf  any  individual  gust 
encountered.  The  data  cannot  be  extended  to  other  flight  conditions 
except  by  assuming  dynamically  similar  conditions. 

Inspection  of  figure  51(a)  indicates  that  the  amplitude  of  the 
rolling  motion  is  much  greater  than  either  pitch  or  yaw.  Similar  results 
are  indicated  in  figure  51(b)  for  the  angular  velocities.  The  data  of 
figure  51  can  be  utilized,  if  dynamically  similar  airplanes  and  a 


58 


NACA  TN  1976 


constant  level  of  atmospheric  turbulence  are  assumed,  to  estimate  the 
fraction  of  the  airplane  life  flown  "before  displacements  or  angular 
velocities  of  the  airplane  in  excess  of  selected  values  are  experienced. 

Relatively  little  data  have  teen  obtained  on  angular  accelerations 
in  rough  air  and  the  data  for  roll  are  discussed  in  the  section  about 
gust  structure.  Twenty-one  values  of  angular  acceleration  in  pitch  were 
available  from  flights  of  the  XB-15  airplane  and  show  that  the  linear 
acceleration  increments  at  the  tail  were  0  to  40  percent  greater  than 
the  acceleration  increment  at  the  center  of  gravity.  The  average 
increase  was  about  25  percent.  If  dynamically  similar  airplanes  are 
assumed,  the  results  can  be  expressed  as  an  average  value  of  the  angular 
acceleration  in  pitch  for  any  airplane  as 


dq  _  20  An 
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At  present  exact  estimations  of  the  angular  acceleration  in  pitch  are 
not  possible. 

The  effect  of  disturbed  motions  on  imposed  loads  is  appreciable 
and  may  affect  the  frequency  distribution  by  as  much  as  50  percent 
(fig.  50).  A  previously  mentioned  analysis  of  piloting  effects  that 
could  be  construed  as  disturbed  motion  indicated  a  variation  of  some 
20  percent  in  the  imposed  loads.  Although  such  variations  cannot  be 
predicted  at  this  time,  current  methods  of  estimating  loads  consider 
this  factor  in  that  the  measured  flight  load  experience  includes  such 
effects. 


RESUME 


A  re'sume  of  the  more  important  findings  under  gust  structure,  air¬ 
plane  reactions,  and  operating  statistics  is  Included  because  of  the 
variety  of  subjects  covered. 


Gust  Structure 

The  gust -gradient  distance  and  gust  size  are  largely  independent  of 
weather,  altitude,  and  airplane  when  expressed  in  terms  of  the  mean 
geometric  chord.  The  gradient  distance  for  a  vertical  gust  has  a 
probable  value  of  10  to  14  chords  for  the  higher  gust  intensities  and  is 
the  same  for  both  spanwise  and  flight  directions. 


NACA  TN  197 6 


59 


The  gust  shapes  vary  widely  hut  triangular  or  sinusoidal  profiles 
in  the  direction  of  flight  seem  good  approximations. 

An  approximation  to  the  spanwise  gust  profile  for  unsymmetrical 
gusts  would  he  a  uniform  distribution  on  which  is  superimposed  a  linear 
variation  of  gust  velocity  with  equal  and  opposite  gust  intensities  at 
each  wing  tip. 

Gust  spacing  for  use  in  estimating  the  number  of  gusts  in  rough 
air  as  a  function  of  gust  intensity  is  about  11  chords  for  velocities 
above  a  threshold  of  0-3  foot  per  second.  The  spacing  is  independent 
of  other  factors.  1 

Lateral,  longitudinal,  and  vertical  gusts  have  essentially  the 
same  characteristics. 

The  frequency  distribution  of  gust  intensities  in  rough  air  is 
independent  of  altitude  and  airplane  when  the  gust  intensity  is  expressed 
in  terms  of  indicated  velocity. 


Airplane  Eeactions 

Aerodynamics. -  Although  the  two-dimensional  unsteady-lift  functions 
yield  the  most  satisfactory  estimates  of  the  load  Increments  due  to 
gusts  for  conventional  airplanes  with  fuselages,  finite-aspect-ratio 
unsteady "lift  functions  appear  pertinent  for  flying  wings.  Tor  unconven¬ 
tional  configurations  such  as  swept  wings,  strip  theory  with  two- 
dimensional  unsteady-lift  functions  appears  adequate. 

The  slope  of  the  lift  curve  for  gust-load  calculation  can  be  best 
estimated  for  engineering  purposes  by  using  an  arbitrary  section  value 
of  6  per  radian  and  simple  aspect -ratio  corrections.  For  swept  wings, 
the  best  estimate  of  lift-curve  slope  is  obtained  by  correcting  the 
value  for  the  equivalent  straight  wing  by  the  cosine  of  the  sweep  angle. 

Available  information  on  compressibility  effects  indicates  that 
the  unsteady-lift  functions  would.be  only  slightly  affected  below  the 
critical  Mach  number  and  would  tend  to  disappear  for  supersonic  speeds. 

At  this  time,  the  use  of  high-speed  wind-tunnel  data  is  suggested  for 
estimating  the  slope  of  the  lift  curve. 

Blgld-body  reactions.-  Analysis  indicates  that  the  mass  parameter 
is  a  significant  measure  of  airplane  reactions  and  that  the  use  of  net 
wing  area  results  in  best  agreement  between  calculation  and  experiment. 
Analyses  neglecting  unsteady-lift  effects  are  not  considered  adequate. 
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No  one  significant  parameter  of  airplane  stability  for  flight  in 
gusty  air  has  been  found.  The  center -of -gravity  position  is  an  important 
factor  for  a  given  airplane  and  arrangement  and  tail  volume  have  some¬ 
what  less  effect.  Comparison  of  experiment  and  calculation  shows  that 
minor  differences  in  the  constants  used  in  calculation  can  vary  the  agree¬ 
ment  from  good  to  poor.  Detailed  calculations  of  airplane  reactions  for 
estimating  loads  are  not  warranted  at  this  time  because  the  accuracy 
required  is  beyond  the  scope  of  present  information. 

Available  data  on  horizontal  and  vertical  tail  loads  indicate  that 
detailed  calculations  of  tail  load  are  not  warranted.  The  lift  on 
either  tail  surface  due  to  a  gust  can  be  estimated  apparently  by 
utilizing  a  true  gust  velocity  and  neglecting  the  alleviating  effects  of 
unsteady  lift  and  airplane  motion  but  considering  the  downwash. 

Elastic -airplane  reactions. -  Analytical  and  experimental  studies 
indicate  that  extensive  simplifications  lead  to  unrealistic  answers  and 
more  exact  solutions  are  required  as  the  wing  mass  becomes  a  greater 
fraction  of  the  total  mass. 

The  effect  of  forward  speed  on  elastic  response  was  negligible, 
but  the  wing-tip  accelerations  increased  with  speed. 

It  does  not  appear  feasible  to  obtain  generalized  solutions  of 
dynamic  response. 


Operating  Statistics 

The  average  miles  to  exceed  limit  load  factor  varies  widely  accord¬ 
ing  to  operator  and  period.  A  most  important  factor  in  setting  the  level 
of  load  appears  to  be  the  probable  speed,  that  is,  the  speed  at  which 
maximum  acceleration  is  most  likely  to  be  experienced,  since  the  gust 
experience  on  all  routes  was  about  the  same.  The  probable -speed  ratio 
(ratio  of  probable  speed  to  the  high  speed  in  level  flight)  is  increasing 
as  new  airplanes  are  introduced  with  a  consequent  reduction  in  the 
average  number  of  flight  miles  required  to  exceed  limit  acceleration 
increment . 

The  flight  miles  to  exceed  the  placard  speed  has  decreased  for  the 
newer  types  of  airplanes. 

The  total  number  of  gusts  encountered  by  an  airplane  during  its 
life  is  a  function  of  the  amount  of  rough  air  flown.  The  amount  of 
rough  air  is  on  the  average  about  10  percent  of  the  total  mileage.  The 
number  of  gusts  encountered  in  rough  air  depends  on  airplane  size  and 
an  average  figure  would  be  500/c  gusts  per  mile  greater  than  0*3  foot 
per  second. 
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The  disturbed  motions  of  an  airplane  in  gusty  air  do  not  correlate 
with  the  maximum  gust  intensities  experienced.  Results  indicate  that 
the  disturbed  motions  may  change  the  loads  in  gusty  air  by  20  to 
50  percent. 


CONCLUDING  REMARKS 


The  available  information  on  gust  loads  does  not  permit  the  deter¬ 
mination  of  applied  loads  under  actual  operating  conditions  on  an  absolute 
basis.  Lack  of  information  on  gust  statistics,  the  behavior  of  airplanes 
in  rough  air,  transient  aerodynamics,  and  operating  conditions  make  such 
an  estimation  of  doubtful  accuracy.  Consequently,  the  procedure  has 
been  to  transfer  loads  data  from  a  reference  airplane  to  new  airplanes 
by  using  available  knowledge  of  aircraft  behavior.  Statistical  data  on 
flight  loads  and  gust  experience  under  actual  operating  conditions  are 
collected  in  order  to  set  the  level  of  loading,  and  "transfer"  coeffi¬ 
cients  based  on  a  knowledge  of  airplane  reaction  are  then  calculated 
and  applied  to  new  airplanes.  No  radical  change  in  operating  conditions 
or  airplane  characteristics  are  assumed  to  be  introduced.  In  addition, 
it  is  implied  that  all • airplanes  are  of  the  same  general  character 
(conventional)  and  that  the  relative  loads  for  single  isolated  gusts  are 
a  measure  of  the  relative  loads  in  a  sequence  of  gusts.  The  important 
features  and  implications  in  this  procedure  are  that  the  same  procedures 
must  be  used  to  evaluate  loads  data  as  are  used  in  loads  calculations. 

In  addition,  the  procedure  presupposes  no  drastic  change  in  airplane 
configuration.  The  use  of  the  same  procedures  of  calculation  tends 
to  eliminate  errors  that  are  due  to  inability  to  predict  airplane 
reactions  on  an  absolute  basis  in  that  errors  in  the  evaluation  of  loads 
data  cancel  those  in  the  calculation  of  loads* 

Extensive  changes  in  configuration  may  be  a  serious  limitation 
since  the  procedure  assumes  in  its  general  application  that  the  effect 
of  pitch  on  gust  loads  is  the  same  for  all  aircraft  and  that  all  air¬ 
craft  respond  to  the  same  gusts  (as  function  of  airplane  size)  to  experi¬ 
ence  significant  loads.  When  new  configurations  such  as  tailless 
airplanes  or  canard  airplanes  are  considered,  the  load  transfer  must 
be  based  on  the  relative  motions  of  the  reference  airplane  and  the  new 
airplane  and  not  on  the  motions  of  the  new  airplane  alone.  A  second 
element,  gust  selection,  has  not  yet  been  resolved  but  it  would  appear 
that  tailless  airplanes,  because  of  small  damping  in  pitch,  and  canard 
airplanes,  because  of  the  adverse  initial  pitching  moment  in  a  gust  for 
swept-wing  airplanes,  might  select  some  other  ranges  of  gust  sizes  than 
'those  selected  by  conventional  aircraft.  So  far,  no  information  is 
available  on  this  problem  and  it  can  only  be  assumed  that  available 
gust -structure  data  are  adequate. 
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In  some  generalized  solutions  of  load  prediction,  a  further  restric¬ 
tion  has  teen  introduced  by  the  sutstitution  of  wing  loading  for  mass 
parameter  as  a  significant  variable.  The  use  of  an  alleviation  factor 
based  on  such  a  substitution  is  still  further  restricted  by  the  assump¬ 
tion  that  all  airplanes  have  the  same  relation  between  wing  loading  and 
mass  parameter.  Information  presented  on  the  load  transfer  coefficient 
indicates  that  the  load  transfer  coefficient  has  appreciable  spread 
when  this  substitution  is  made  and  the  spread  seems  to  be  a  function  of 
airplane  type  and  use  since  the  data  can  be  grouped  according  to  land 
transports,  flying  boats,  and  personal  airplanes.  If  wing  loading  is 
used,  therefore,  different  relations  should  be  considered  according  to 
airplane  class  and  use. 

In  cases  where  the  dynamic  response  may  be  a  factor,  it  should  be 
considered  as  an  amplifying  factor  for  the  incremental  loads  whether 
for  studies  of  large  loads  or  repeated  loads.  Eynamlc  response  should 
be  considered  on  a  relative  basis  since  the  uncertainties  in  calcula¬ 
tions  due  to  lack  of  knowledge  and  practical  simplifications  have  not 
been  resolved. 

Of  the  many  related  problems,  that  of  gust-alleviating  systems  is 
of  great  importance.  Such  systems  are  of  particular  interest  in  reducing 
loads  and  improving  riding  comfort  for  high-speed  aircraft.  The  two 
problems  may  not  have  compatible  solutions  since  'load  reduction  implies 
the  reduction  of  stresses  in  all  critical  members  without  impairing 
other  qualities  of  the  airplane  while  the  improvement  in  riding  comfort 
implies  the  reduction  of  airplane  motions  and  accelerations.  Although 
the  basic  problem  of  reducing  the  load  due  to  a  single  gust  is  amenable 
to  analysis,  the  related  problems  of  maintaining  other  airplane  qualities 
and  of  determining  alleviation  in  rough  air  introduce  difficulties. 

Three  fundamental  systems  of  gust  alleviation  are  available  and 
consist  of  changing  airplane  characteristics,  aeroelastic  systems,  and 
aeromechanical  systems.  The  first  method  consists  of  emphasizing  the 
pertinent  characteristics  of  the  airplane,  such  as  increasing  the  wing 
loading  to  reduce  accelerations,  decreasing  the  slope  of  the  lift  curve 
(by  use  of  vents  or  sweep  and  by  changing  the  aspect  ratio) ,  adjusting 
the  stability  characteristics  to  obtain  favorable  response  in  rough  air, 
or  using  free-floating  flaps.  Aeroelastic  systems  generally  involve 
torsional  deflections  of  the  wing  due  to  imposed  loads.  Such  systems 
can  be  designed  by  incorporating  slant  hinges  or  flaps  operated  by  wing 
deflections  or  by  the  adjustment  of  the  torsional  and  bending  rigidity. 
Aeroelastic  systems  have  the  advantage  of  being  inherent  in  the  construc¬ 
tion  of  the  airplane  and  the  disadvantage  that  they  form  an  additional 
source  of  elastic  phenomena  such  as  flutter  or  dynamic  response  due  to 
a  series  of  gusts.  Aeromechanical  systems  cover  combinations  of  servo¬ 
mechanisms,  detectors,  and  aerodynamic  controls.  Some  detector  systems 
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considered  are  the  use  of  an  accelerometer  or  strain  gage  to  measure 
accelerations  or  load  and  the  use  of  a  feeler  vane  or  similar  system  to 
detect  gust -velocity  or  angle-of-attack  changes.  The  servo-systems  are 
assumed  to  operate  the  conventional  controls  of  the  airplane  or  to 
operate  wing  flaps  or  spoilers.  The  difficulties  are  those  of  possible 
failures  of  complicated  mechanisms,  the  question  of  oscillation  of  the 
system,  the  response  of  repeated  gusts,  and,  finally,  the  need  as  air¬ 
plane  speed  increases  of  obtaining  a  servo-system  that  has  a  high  fre¬ 
quency  response  (imposed  frequencies  up  to  10  cps) .  The  advantages 
are  flexibility  and  utilization  of  many  of  the  available  mechanisms. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base,  7a.,  August  5,  19h9 
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APPENDIX  A 

COOPERATING  AIRLINES  AND  AGENCIES 

The  following  airlines  and  agencies  have  cooperated  extensively 
with  the  NACA  in  obtaining  much  of  the  data  used  in  the  preparation  of 
this  paper: 

American  Airlines,  Inc. 

Eastern  Air  Lines,  Inc. 

Pan  American  World  Airways  System 
Trans  World  Airline,  Inc. 

United  Air  Lines,  Inc . 

Department  of  Commerce 
U.  S.  Weather  Bureau 
Civil  Aeronautics  Administration 
U.  S.  Air  Force 

Air  Materiel  Command 
Air  Weather  Service 

U.  S.  Navy,  Bureau  of  Aeronautics,  Structures  Unit 


O  I 
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APPENDIX  B 


SYMBOLS 


A 

B 

L 

C 


aspect  ratio  (b^/S) 

arbitrary  constant 

wing  span,  feet 

arbitrary  constant 

mean  geometric  chord  of  wing,  feet 

mean  aerodynamic  chord.,  feet 

lift  coefficient 

variation  of  lift  coefficient  when  penetrating 
a  sharp-edge  gust,  expressed,  as  a  fraction 
of  final  lift  (Kussner's  function,  refer¬ 
ence  2) 

variation  of  lift  coefficient  following  a  sudden 
change  in  angle  of  attack,  expressed  as  a 
function  of  final  lift  (Wagner's  function, 
reference  2) 

pitching -moment  coefficient  about  center  of 
gravity 

differential  operator 


4S>g/dCL 

f  ■' 


g 


static  margin 

frequency  of  occurrence 

wing  frequency,  cycles  per  second 

acceleration  due  to  gravity,  feet  per  second  per 
second 
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H 

K 

L 

1 

M 

n 

An 

Ans 


Anr 

P 

q. 

s 

s 

t 

u 

ue 


gust -gradient  distance  (distance  in  which  the 
vertical  velocity  of  a  gust  rises  linearly 
from,  zero  to  its  maximum,  value) ,  feet  or  chords 

gust  alleviation  factor 

lift,  pounds 

distance  from  trailing  edge  of  wing  to  leading 
edge  of  tail,  chords 

Mach  number 

normal  acceleration,  g  units 

acceleration  increment  due  to  a  gust,  g  unit 8 
(n  -  1  for  most  cases) 

acceleration  increment,  g  units 

^computed  according  to  formula 

acceleration  increment  on  rigid  airplane,  g  units 
probability 

dynamic  pressure,  pounds  per  square  foot 
wing  area,  square  feet 

distance  airplane  has  penetrated  gust  at  time  t, 
chords 


time  from  start  of  gust  penetration,  seconds 
(cs/V) 

gust  velocity  at  peak  of  gust,  feet  per  second 
effective  gust  velocity,  feet  per  second 

( computed  from  recorded  acceleration  increment 
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u 

V 

Ya 


vmax 

Vnev 

V„ 


W 
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vertical  gust  velocity  at  axiy  point,  feet  per 
second 

true  airspeed,  feet  per  second 

equivalent  airspeed,  feet  per  second 

probable  speed,  miles  per  hour 

design  level-fligjit  speed,  miles  per  hour 

maximum  indicated  airspeed,  miles  per  hour 

placard  do-not -exceed  speed,  miles  per  hour 

airspeed  at  maximum  acceleration  increment, 
miles  per  hour 

weight  of  airplane,  pounds 


z 


vertical  displacement  of  airplane,  feet 


dp/dt 


angular  acceleration  in  roll,  radians  per  second 
per  second 


dq/dt 


angular  acceleration  in  pitch,  radians  per  second 
per  second 


de/da 


dovnwash  factor 


angle  of  attack,  degrees 

angle  of  gust  to  direction  of  flight  of  airplane, 
degrees 

incremental  value 

absolute  displacement  of  equivalent  fuselage,  feet 

absolute  displacement  of  equivalent  wing,  feet 

deflection  of  equivalent  wing  under  conditions 

corresponding  to  normal  static  design  procedure, 
feet 

maximum  value  of  5^  -  5p ,  feet 


0 


pitch  angle,  degrees 
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gust  spacing  (distance  "between  gust  peaks),  chords 


mass  parameter  (|ig  is  twice  the  relative  density 
(reference  2k)  divided  "by  the  lift-curve  slope) 


pgSc  ZT 

CLCX  1 


Subscripts: 


density  of  air,  slugs  per  cubic  foot 

density  of  air  at  sea  level,  slugs  per  cubic  foot 


average 

probable;  denotes  the  most  probable  value  of  a 
quantity  as  determined  from  experimental  data 

denotes  value  computed  from  local  measurements 

center  of  gravity 

wing 

for  extended  equations  denotes  point  at 
which  solution  applies 


maximum 


Subscripts  used  with  An: 


acceleration  increment  due  to  a  gust 

acceleration  increment  due  to  vertical  motion  of 
airplane 

acceleration  increment  due  to  pitch  of  airplane 

acceleration  increment  due  to  angular  velocity  of 
airplane 

acceleration  increment  due  to  downwash 


acceleration  increment  on  wing 
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s 

T 


Examples: 


acceleration  increment  on  stabilizing  surface 
total  acceleration  increment 


acceleration  increment  on  wing  due  to  gust 
total  acceleration  increment  due  to  tail 
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TABLE  V 

FREQUENCY  OF  OCCURRENCE  OF  GIVEN  VALUES  OF  ANGULAR  ACCELERATION  » 

IN  ROLL  WITH  ASSOCIATED  VALUES  OF  NORMAL 
ACCELERATION  INCREMENT 

(XC-35  airplane  data] 
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TABLE  VII 


CHARACTERISTICS  OF  AIRPLANE  MOXiEL 


jjamlnar-flow  wingj 


Weight,  It . 11. 84 

Wing  area,  sq.  ft . 5.44 

Wing  loading,  lt/sq  ft . 2.l8 

Span,  ft  . .  7 

Mean  aerodynamic  chord,  ft . O.829 

Aspect  ratio  .  .  9 

Taper  ratio . O.382 

Center  of  gravity,  percent  M.A.C . 30-5 

Wing  section . 65-3 “4l8,  a  =  1-0 

Lift -curve  slope,  per  radian 

Smooth  wing  (steady  flow)  . 5-01 

Roughened  wing  (steady  flow) . 4.l8 

Theoretical . 5*^1 


COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED 
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CHARACTER  1ST  ICS  OF  GUST  -TUNNEL  MODELS 
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TABLE  XII 


MAXIMUM  TOTAL  ACCELERATION  INCREMENTS  ON  CONVENTIONAL  AIRPLANE 
[All  increments  given  in  g  unitsj 


(a)  Center  of  gravity,  15  percent  M.A.C. 
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TABLE  XII  -  Continued 


MAXIMUM  TOTAL  ACCELERATION  INCREMENTS  -  Continued 


("b)  Center  of  gravity,  25  percent  M.A.C. 
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TABLE  XII  -  Concluded. 


MAXIMUM  TOTAL  ACCELERATION  INCREMENTS  -  Concluded 


(c)  Center  of  gravity,  35  percent  M.A.C. 


Tall  length  Tall  area  Aru  Aru  An,  Ail,  An™  Aru  Aru  Au  An,  Anfi  An  Ad™  An™ 
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TABLE  XVI 


TEST  MODEL  CHARACTERISTICS 


[Nonrigid  airplanej 


Weight,  lh  . 

Wing  area,  sq  ft . 

Mean  geometric  chord,  ft  ..... 

Span,  ft . . 

Slope  of  lift  curve,  per  radian  •  . 

Forward  velocity,  fps  . 

Gust  velocity,  fps  ........ 

p 

Pitching  moment  of  Inertia,  slug-ft 
Radius  of  gyration  of  wing,  ft  .  . 
Weight  of  wing,  lh  . 


.  I.832 

•  1.183 

•  0.394 

.  3-0 

.  4.73 

.  6l.O 

6.0 
0.00782 
.  0.663 


0.293 
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denotes  probable  speed  of  flight. 
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Figure  1.—  Variation  of  normal  acceleration  with  angle  of  gust  for  four  values  of  wing  loading 

Y  =  200  miles  per  hour;  U  =  15  feet  per  second. 
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/vase  parameter, pig 

Figure  4.—  Acceleration  ratio  as  function  of  the  mass  parameter.  Above 
dash  line,  finite  gradient  solution  is  not  valid. 
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Figure  5.—  Variation  of  alleviation  factor  with  wing  loading. 

From  reference  1. 
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Altitude,  ft 
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;ure  8.—  Probability  distribution  defining  Figure  9*—  Probability  that  an  effective  gus- 
limits  of  gust  distributions.  From  velocity  at  various  altitudes  will  exceed 
reference  9-  a  given  value.  XC— 35  airplane.  From 

reference  10. 
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Figure  11.—  Gust  velocity  U  as  a  function  of  gust-gradient  distance 
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(b)  Gust— gradient  distance  in  spans. 

Figure  12.—  Continued. 
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Figure  13 Belation  between  the  probable  and  average  gust-gradient  distance. 
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Figure  15.—  Span-wise  gust  distributions  and  frequency  of  occurrence  f 

for  each  type. 
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>ur®  17 •—  Comparison  of  the  effective  gust  velocity  obtained  from  acceleration  records  vith  the 
average  effective  gust  velocity  along  the  span  of  the  XC— 35  airplane.  (Data  chosen  at  random 
from  25  percent  of  the  total  number  of  flights.)  From  reference  12. 
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jure  20.-  Comparison  of  the  maximum  vertical  Figure  21.-  Variation  of  vertical  gust  velocity 

gust  velocity  with  the  maximum  horizontal  gust  with  horizontal  gust  velocity  for  equal 

velocity  for  separate  traverses  in  rough  air.  frequencies  of  occurrence. 

Data  for  XC— 35  airplane. 
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(a)  Angle— of -at  tack  change. 


(b)  Lift-coefficient  change. 

I’igure  22.—  Diagrammatic  illustration  of  the  superposition  (Duhamel) 
integral  for  evaluation  of  lift  following  arbitrary  variation  in 
angle  of  attack. 
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Figure  25.—  Unsteady— lift  function  for  an  airfoil  traversing  a  sharp  gust 
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Figure  29.—  History  of  acceleration  for  an  airplane  with  a  45°  sweptback  wing  in  a  sharp— edge  gust 
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Figure  31 •—  Test  model. 
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Figure  32.—  Experimental  acceleration  ratio  for  Boeing  B— 247  airplane  as 

a  function  of  gust-gradient  distance. 


Figure  33-—  Effect  of  tail  length  on  the  -vertical  Telocity  at  the  tail  surface  following  a  sudden 

increase  in  wing  circulation.  From  reference  31. 
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(a)  H  =  0.  (h)  H  =  8.  (c)  H  =  16. 

Figure  34.—  The  ratio  of  total  wing  load  to  the  wing  load  with  no  pitch  as 
a  function  of  static  margin  for  a  conventional  airplane. 
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Figure  35*—  Comparison  of  calculated  and  experimental  acceleration 
increments  for  conventional  airplane.  All  points  except  where 
noted  have  center  of  gravity  at  25  percent  mean  geometric  chord. 
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Figure  36.—  Variation  of  relative  lift  due  to  pitch  with  masB  parameter 
for  a  gust-gradient  distance  of  8  chords. 


(a)  E  =  0.  (b)  H  =  8.  (c)  H  =  16. 


Figure  37*—  Tail  load  increment  as  a  fraction  of  tail  load  due  to  gust 

only  for  a  conventional  airplane. 
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(a)  Tailless  airplane. 


(t> )  Canard  airplane . 

Figure  38.-  Comparison  of  calculated  and  experimental  acceleration 
increments  for  unconventional  airplanes. 
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Figure  39 •—  Relative  response  of  representative  airplanes  as  a  function 
of  ving  loading.  (Acceleration  based  on  W/S  divided  by  acceleration 
based  on  mass  parameter.) 


Figure  40.—  Relative  response  of  an  airplane  encountering  sinusoidal  and 
triangular  gusts  as  a  function  of  mass  parameter. 
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Figure  ■4-3.—  Comparison  of  calculated  and  experimental  values  of  ratio  of  maximum  ving—t  ip— deflection 
increment  to  maximum  acceleration  increment  for  change  in  gradient  distance. 
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Figure  47.—  Effect  of  change  of  ving  frequency  due  to  addition  of  mass 

Model  D,  conditions  1  and  2. 
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Figure  48.—  Fraction  of  flight  time  spent  at  speeds  equal  to  or  greater  than  selected  values 

during  airline  operations. 
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Hire  49.—  Average  time  required  to  exceed  a  Figure  50.  Influence  of  disturbed  motion  of 
given  value  of  airspeed.  Airplane  E.  XC“35  airplane  on  apparent  gust  frequency 

distribution  in  very  rough  air. 
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(b)  Relative  frequency  of  given  values  of  angular  velocity. 

Figure  51*—  Relative  frequency  with  which  selected  values  of  angular 
motions  will  he  exceeded  in  rough  air. 
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